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1 Introduction 

The Allen-Cahn equation 

— Au = u — u 3 , |it| < 1 (1) 

has been studied for several decades and is an important nonlinear PDE due to 
the fact that it lies at the interface of several different mathematical fields. The 
famous De Giorgi conjecture states that any entire solution to (JT]) in I" which 
is monotone in one direction should be one dimensional, at least for n < 8. 
The conjecture was proved in dimension n = 2 by Ghoussoub-Gui (USD and 
dimension 3 by Ambrosio-Cabre (in), and in dimensions 4 < n < 8 by Savin 
dm under an additional assumption. For n > 9, counter-examples have been 
constructed by del Pino-Kowalczyk-Wei (El)- Note that monotone solutions 
are indeed minimizers with respect to local perturbations (HD- 

A natural extension of De Giorgi’s conjecture is to classify stable or finite 
Morse index solutions. Regarding stable solutions, it has been shown (in. m) 
that stable solutions in R 2 are necessarily one-dimensional, while Pacard-Wei 
(p4]l constructed stable solutions in M 8 which are not one-dinrensional. 
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The study of finite Morse index solutions is much more involved. In R 2 , 
we have now a rather complete picture of Morse index one solutions. They are 
so-called four-end solutions, parametrized by the angle between the two lines. 
The cross solution, constructed by Dang-Fife-Peletier ([8]), represents a four- 
end solution with angle j, while the almost parallel line solution, constructed 
by del Pino-Kowalczyk-Pacard-Wei (EH). represents a four-end solution with 
angle close to ^ or 0. The existence of four-end solutions with any angle between 
0 and ^ was proved by two methods: the first through the moduli space theory 
by Kowalczyk-Liu-Pacard (isi). and the second approach by the mountain-pass 
variational method by us (Gui-Liu-Wei [T5]). It was also shown that the four- 
end solutions have Morse index one ([IS])- On the other hand, monotone and 
symmetric properties of a general four end solution have been obtained in [14] , 
where more general finite morse index solutions in M 2 have also been studied 
under an extra energy condition or a condition on the asymptotical structure of 
nodal curves. 

In this paper we are interested in the structure of two-end solutions of m 
in R 3 . It turns out that without the monotone condition, there are actually 
a lot of solutions. One simple example is the so called saddle solution whose 
nodal sets are precisely the xoy,yoz and xoz planes (Alessio-Montecchiari [3]). 
del Pino-Kowalczyk-Wei ([T2]) proved that for each non-degenerate minimal 
surfaces with finite total curvature, one could find a solution to m whose nodal 
sets are close to a rescaled version of this minimal surface. In particular, there 
are axially solutions whose nodal sets are close to catenoids with very large 
waist. Axially symmetric solutions with multiple interfaces which are governed 
by the Jocobi-Toda system are constructed in Agudelo-del Pino-Wei [2]. 

In spite of all these developments, some important questions for m in R 3 
remain unanswered, even for axially symmetric solutions. In this paper, we will 
study those axially symmetric solutions which are additionally even with respect 
to the xoy plane. In terms of the cylindrical coordinate (?’, z) , they satisfy 

f u zz + u rr + r~ 1 u r + u — u 3 = 0, r £ [0, +oo), z £ M, , . 

I u (r, z) = u (r, —z ), (0, 2 ) = 0. ' 

Let H (x) = tanh -^= be the one dimensional heteroclinic solution: 

-H" = H-H 3 ,H (0) = 0, H (± 00 ) = ±1. 

The solutions we are interested will have H as its asymptotic profile. We say 
that a solution u of m has growth rate k if it has the following asymptotic 
behavior: 


ll M (r, ■) — H (• — k\nr + c)|| ioo(0|+oo) 0, as r -s> + 00 , (3) 

for certain c £ R. The existence results obtained in [5] and [12] based on 
Lyapunov-Schmidt reduction arguments tell us that there are solutions whose 
growth rate is in the interval (V2, \[2 + <5) and (<5 -1 , + 00 ) , where S is a very 
small constant. A natural question is, whether or not there are solutions with 
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growth rate in the range [-s/2 -f- 5,5 x ] . In this paper we answer this question 
affirmatively and our main result is the following 

Theorem 1 For each k £ (s/2, +oo) , there exists a solution to @ which has 
growth rate k. 

As we will see later, these solutions indeed are monotone in the following 
sense: 

u z > 0 for z > 0; u r < 0, for r > 0. (4) 

Outside a large ball, the nodal set of the solutions given by Theorem [Q has 
two components, each component is asymptotic to a catenoidal end and around 
each end, the solution look likes the one dimensional heteroclinic solution. Bor¬ 
rowing a terminology from minimal surface theory, a solution satisfying ED and 
ED will be called a two-end solution. We emphasize that here by definition the 
two-end solutions are axially symmetric. Comparing with the corresponding 
definition of two-end minimal surfaces, it seems that a more general definition 
of two-end solution should also involve those solutions which are not axially 
symmetric and only assume that their nodal set is asymptotic to two catenoidal 
ends. In the minimal surface theory, a classical result proved by R. Schoen ([28]) 
is that a minimal surface with two catenoidal ends is a catenoid. We expect that 
the analogous result for Allen-Cahn equation should hold. A major difficulty 
we encounter is to show that a solution whose nodal set is asymptotic to two 
catenoidal ends is axially symmetric. 

Observe that for each k £ (0, + 00 ), there exists a catenoid with growth rate 
k. Taking into account the relation between minimal surface theory and Allen- 
Cahn equation, at first glance, one may think that for each k £ (0, + 00 ), there 
should be a two-end solution of Allen-Cahn equation. But this turns out to be 
false. In fact, we have 

Theorem 2 There does not exist two-end solution with growth rate k £ (0, ^j~\- 

Indeed, one expects that for k £ (^,y/2], two-end solution with growth 
rate k also should not exist, while for each k £ (x/2, + 00 ) , there should be a 
unique two-end solution with growth rate k. We remark that the lower bound 
x/2 is somehow related to the deep facts that the two ends of the solutions to 
the Allen-Cahn equation actually “interact” with each other and in this regime 
one naturally encounters the so called Toda system, as we will see later in the 
analysis. This constitutes a major difference with the theory of minimal surfaces. 
Roughly speaking, the Allen-Cahn equation interplays between the theory of 
minimal surfaces and the theory of Toda system, which is a classical integrable 
system. In the minimal surface theory, the catenoids are basic blocks for the 
construction of other minimal surfaces or constant mean curvature surfaces (see 
for example m, m , m and the references therein). It is therefore natural 
and interesting to ask what role two-end solutions of the Allen-Cahn equation 
play in constructing other solutions. We also remark that the two-end solutions 
given by Theorem[l]are all unstable, and the stable solution conjecture says that 
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the only bounded stable solution of the Allen-Cahn equation in R 3 should be 
one dimensional. It is also expected that the Morse index of two-end solutions 
should be equal to one, again, we don’t have a proof of this statement, although 
it is known that the two-end solutions constructed in [2] and m have Morse 
index one. 

The results in Theorem |T| could be regarded as a generalization of the cor¬ 
responding results for four-end solutions in R 2 ([18], BSD- To explain this, let 
us say a few words about the multiple-end solutions in R 2 . By definition, a 
2fc-end solution of m in R 2 is a solution whose nodal set outside a large ball is 
asymptotic to 2k half straight lines at infinity. It is known (del Pino-Kowalczyk- 
Pacard 0) that the set of 2fc-end solution in R 2 has a structure of real analytic 
variety of formal dimension 2k. Some examples of solutions near the boundary 
of this “moduli space” have been constructed (del Pino-Kowalczyk-Pacard-Wei 
m, Kowalczyk-Liu-Pacard-Wei [22]) ■ As we mentioned above, for k = 2, it 
is proved in Kowalczyk-Liu-Pacard [T9: that the moduli space of four-end solu¬ 
tions, modulo rigid motion, is diffeomorpic to the open interval (0,1). It is also 
proved there that for each 9 € (0, f) , there exists a four-end solution ug which 
is even with respect to both x and y axis and the asymptotic line of the nodal 
set of ug in the first quadrant makes angle 9 with the x axis. Now we observe 
that if we reflect the solution to (0 across the z axis, then we get a solution 
defined for all (r, z) 6 R 2 and even in both variables. Moreover, the nodal set of 
this solution outside a large ball also has four components. Hence the two-end 
solutions in R 3 are in certain sense analogy of the four-end solutions in R 2 . The 
equations they satisfied are different from each other only by the term r -1 u r 
which makes the problem inhomogeneous. In fact in R 2 a major fact used is 
that there are two linearly independent kernels u x and u y . The construction of 
four-end solution is done by first proving that all four-end solutions are nonde¬ 
generate. Here we face the problem of degeneracy. We overcome this difficulty 
by applying global bifurcation theory developed by Buffoni, Dancer and Toland 
for real analytical variety of formal dimension 1. 

Before proceeding to proof of our main results, let us outline the main ideas 
of the proof, since this also gives a description of the set of solutions in Theorem 
[TJ The main steps of the proof are similar to the corresponding analysis for four- 
end solution in dimension two, performed in m and M- Our basic strategy 
is to analyze the structure of the set M of axially symmetric two-end solutions 
with growth rate larger than y/2. In section [21 we prove that the solutions in 
M with some additional natural constraints are compact in certain sense. In 
section [3l we show that M has the structure of a real analytic variety of formal 
dimension 1. In section [4l we analyze those solutions near the “boundary” of the 
moduli space M and prove that they are unique in certain sense. Combining 
these results, we conclude the proof by applying a structure theorem for real 
analytic varieties. We remark that there are two main differences between the 
R 2 and R 3 case. Firstly, the proof of compactness in 3D case is much more 
delicate than the 2D case and detailed asymptotic analysis is needed. Secondly, 
as we mentioned, the four-end solutions in 2D are all non-degenerate (BZI> BSI). 
but we don’t know whether it is true for two-end solutions in 3D. 
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2 Compactness of two-end solutions 

The compactness of moduli spaces of minimal surface plays an important role 
in the minimal surface theory, for example, it is an important step towards the 
classification of certain type of minimal surfaces. We refer to Perez-Ros [55] and 
references therein for more details on this and other related subjects. In this 
section, we shall investigate the compactness property for two-end solutions of 
the Allen-Cahn equation. 

Throughout the paper we denote by E the set [0, +oo) x M and by E + the 
set [0, Too) x [0, Too). Let {n ra } be a sequence of two-end solutions which has 
growth rate k n > y/2. (Note that for k £ (^,-\/2], we don’t know whether or 
not there is a two-end solution with growth rate k.) Then by definition, 

||w„ (r, •) - H (• - k n Inr- - c n )|| iao(0 +oo) —> 0, as r ->• Too. (5) 

We will show that if the distance of the nodal set of u n to the origin is uniformly 
bounded with respect to n, then up to a subsequence, {u n } converges strongly 
to a two-end solution u x - Here converging strongly means that u n —> u a 0 in 
Cioc (IE) and there exist constants koo, Coo such that k n —> k x , c n —> c x , 

11 Uoo (r,-) - H (• - koo lnr - Coo)|| L oo (0i+00 ) 0, as r ->■ Too, 

and 

u n — H {z — k n In r - c n ) -)■ Uoo - IL (z - fcoo In r - Coo) 

in L°° (E). 

Under the assumption that a solution u has the asymptotic behavior ([5]) , 
one could actually prove that u has certain monotonicity property. 

Lemma 3 Suppose u is a two-end solution with growth rate k > \[2. Then 

d r u < 0 for r > 0; d z u > 0 for z > 0. (6) 

The proof of this result is based on the moving plane method and its proof 
will be given in the appendix. By Lemma [3l the two-end solutions we are 
analyzing always satisfy © . We will denote the nodal set of a solution u in the 
upper r-z plane by 

A f u := {p G E + , u (p) = 0} . 

Due to the monotonicity property, the set J\f u D <9E + contains a unique point, 
call it V u - 

In the rest of the paper, we use C and a to denote universal constants which 
may vary from step to step. The main result of this section is the following 

Proposition 4 Let u n he a sequence of two-end solutions with growth rate 
larger than \[2. Assume 

\V U J <C. 

Then there exists a two-end solution Uoo such that up to a subsequence {u„} 
converges strongly to Uoo. 
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The rest of this section is devoted to the proof of Proposition [I] 

By Lemma [3j for each solution u n , J\f Un will be the graph of a function 

z = fn{r) ,r € [t„,+oo). 

Here t n satisfies V Un = (t ra , f n (t„)). In particular, if V Un is on the z axis, then 
= 0 . 

Lemma 5 Under the assumption of Proposition^ we have 

lim f n (r) = +oo, 

r—>•+oo 

uniformly in n. 

Proof. We argue by contradiction. If this was not true, there will exist a 
constant Co such that for any l £ N, one could find a solution u ni satisfying 

fm (r) < C 0 , for r £ {t n „l). 

Since \u n \ < 1, the sequence {u ni } will converge in Cf oc { E) to a nontrivial 
solution W whose nodal set is contained in the strip \z\ < Co- Moreover, by 
the monotonicity of u n , W r < 0 for r > 0. The limit w (z) := lim r -s.+oo W (r, z) 
then exists and is a solution of the Allen-Cahn equation in dimension 1 : 

— w" = w — w 3 , in R. (7) 

By the symmetric property of W, w (z) = w (—z). We also have w ( z) 1, as 

z —> Too. But m does not have a solution with nodal set containing only in a 
bounded set {\z\ < Co}. This is a contradiction and the proof is finished. ■ 
By the monotonicity of u n , J\f Un is also the graph of a function over the z 
axis: 

N Un = {(r,z) :r = g n (z)} . 

Similar arguments with slight modification as that of dSJ) imply that 

lim g n (z) Too, 

z—*-Y oo 

also uniformly in n. 

To obtain more information on the functions /„, we should use the balancing 
formula (El)- Let A' = (0,0,1) be the constant vector field on M 3 . If u = 
u ( x , y , z) is a solution to the Allen-Cahn equation, then one could check that 

div | ( ^ |Vu| 2 T F {u)^j X - (Vti • X) Vuj = 0. 

Here F (u) = j ( u 2 — l) 2 . Therefore for each regular domain f2, we have the 
following balancing formula: 

/ {(I I v «| 2 + ^ («)) JT - ( y « - JT) v „} - ^ = °. (8) 

Here v is the outward unit normal vector of the boundary dQ. 

We would like to use JH} to control the slope of the function /„. 
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Lemma 6 For each e > 0, there exists t e > 0 such that for all n £ N, 


f'n O 2 ) ~ fn Oh) < e for t e < n < r 2 . 

Proof. We first show that for each S > 0 and f > 0, there exists r* > f, such 
that 

fn ( r n ) < <5, n G N, 

for some r n € ( r,r *). Indeed, if this was not true, then for any f > f, there 
exists n, such that 

fn ( r) >S,r& (f,f). (9) 

Let r < r\ < r 2 < f. Consider the region C K 3 given by 

{(x, y, z) : z > 0, L x (r) < z < L 2 (r)} , 


where 


Li (r) = f n (n) - ' (r - ri), 

fn Oh) 

^2 (r) = /„ (r 2 ) - ) , (r - r 2 ). 

fn M 


By the balancing formula, 


J | Q |Vu| 2 + F (u)J X - (Vu • X) Vu j • xds = 0. 

Using m and the fact that u exponentially decays to 1 away from the interface, 
we deduce that as ri —> +oo, 


J 5r2n{2=o} 
On the other hand, 


l |Vzt| 2 + F (u) ) X - (Vu • X) Vu } ■ vds 0. 


/ d£ir\{z=Li(r)} 

rifn (n) 


i |Vu| 2 + F (u) ) X - (Vu • X) Vzi • vds 


-ci 


\A + /n(r l) 2 


where ci is a fixed constant, and 


t d£LC\{z=L2{r)} 

r 2fn {r 2 ) 


i |Vzt| 2 + F (u) ) X - (Vu • X) Vw !> • vds 


Cl 


\Jl + fni^f 


7 








Combining these estimates, we get a contradiction if r 2 is large enough. 

Fix an e > 0, if the conclusion of the lemma was not true, then for any l > 0, 
one could find l < r± < r 2 such that 


f'n M - f'n (n) = e 


( 10 ) 


and 

1 3 

2 £ < I f'n (01 < 2 £,r £ ( ri ’ r2 ) ' 

In this case, similarly as above, we could estimate 


J dQn{z=o} 
= o(i)n. 


=- |Vu| 2 + F(u))X- (V« • X) Xu !> • vds 


Therefore the balancing formula tells us that 


f'n (f 2 ) f"2 ~ f'n (fl) =0(l)n. 


This contradicts with m- ■ 

Lemma [6] tells us that the slope of the nodal line could not increase much as 
r increases. This in particular implies the following 

Lemma 7 Under the assumption of Proposition^ we have 


lim f n (r) = 0, 

r—H- oo 


uniformly in n. 

Proof. We argue by contradiction and assume that there exists 6 > 0 such that 
for each Co > 0, one could find to > Co and n satisfying 

fn > S. 

Then one could find t \, f 2 > such that 

fh(r) e 

fn(tl) = ^,fn(t 2 ) = S. 

This contradicts with Lemma [G] ■ 

Since | « ra | < 1, one could assume without loss of generality that u n —> Uoo 
in Cf oc (E) for a solution Note that u could not be identically 0. By the 
nronotonicity property of u n , is also monotone, that is, satisfies (JJ). We 
need to prove u n —> in the strong sense, which in particular implies that u ac 

is a solution to ([2]) whose nodal line is asymptotic to a log curve at infinity and 
therefore a two-end solution. 
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The previous lemmas yield certain information on the nodal curve A f Un , but 
we still don’t have precise description of their shape. On the other hand, we 
know that in the region where r is large, locally around the nodal line, the 
solution resembles the heteroclinic solution. Our main idea to prove the strong 
convergence of {u n } is to get a precise decay estimate of the solution u n to H 
along their nodal lines. For notational simplicity, sometimes we will drop the 
subscript n for the function u n , /„. 

To proceed, let us introduce the Fermi coordinate around the smooth curve 
z = f (■) . This coordinate, denoted by (r i, Z\) , is defined through the relation 

(r,z) = (ri,/(n)) + 2 in, 


where n is the unit normal vector at the point (ri, f (ri)) of the curve z = f (■), 
upward pointed. Explicitly, 


l Z = f + Z '7^F' 


(ii) 


Here the function / is evaluated at the point r i. The map X : (ri,zi) —> 
(r, z) from the Fermi coordinate to the original coordinate system will be a 
diffeomorphism in certain tubular neighborhood of the graph of /. 

From (fill) , we get 


f d r ri 

d z ri 

\ - 1 

f d ri r 

d Zl r \ 

V d r zi 

d z zi 

1 ~ \ 

v d ri z 

d Zl z ) 


-l 



Vi+F- 

y/l+f ' 2 


(1 + f r2 )B 

f _ 

V 1 +f 2 V^r- 


f 

(1 + f ,2 )B 

1 


( 12 ) 


Here 

„ Z\f" 

B = 1--- r . 

(i + r 2 r 

Let us denote by A( r>z ) = + d\ and A( rijZl ) = d’^ +9 Zl . These two operators 

are related by 


where 


r,z) A(ri,*i) + d ri 

1 d Zl A 1 d ri A 
+ 2— d «-^ dr " 


A = 1 + f' 2 


f 


f" 2 

(i + /' 2 ) 2 ' 


(13) 


(14) 
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This follows from direct computation, see m for more details. From formula 
(fl3l) . one sees that if is small, then A( r>2 ) is a small perturbation 

of One should keep in mind that in (1131) . there are terms (in d ri A) 

involving the third derivatives of / and these terms should be dealt with very 
carefully. 

For any function 0(r, z), the pullback of O by the diffeomorphism X will 
be defined as 

X*0 (n, zx) := 0 o X (n, zi). 

Occasionally, X*0 will also be denoted by 0*. Keep in mind that X*0 is only 
defined in the region where the Fermi coordinate is well defined. Conversely, for 
any function 0 (n, Zi), we set (A -1 ) 0 (r, z) = 0 o X~ l (r, z). 

To prove the compactness, we need to know the precise asymptotic behavior 
of u n . We shall define suitable approximate solutions and analyze the error 
between the approximate and true solutions. It turns out that to get a good 
approximate solution, it will be convenient to use the Fermi coordinate. Clearly 
there is a technical issue concerning the Fermi coordinate. Namely, the Fermi 
coordiante is in general not well defined in the whole r-z plane. Note that 
for r\ large, by Lemma [7] /' is small. We also know that u is close to the 
one dimensional heteroclinic solution locally around the nodal line at far away. 
Hence intuitively, for r large, f" ..., should also be small. Indeed, we have 

Lemma 8 isr tends to infinity, \ f" (r)|, | (7’) | , |/^ 4 ^ (r)| tend to zero. 

Proof. Let r 0 be a fixed large constant. Around the point (r 0 ,/(?'o)), let us 
consider the line l\ 

z-f{r 0 ) = f (r 0 )(r-r 0 ). 

Let I 2 be the line orthogonal to it and passing through (r 0 ,/(ro )). Use these 
two lines to build an orthogonal coordinate system, denoted by (s,f), where 
the s coordinate corresponding to line l\. We know that u is close to H ( t ). By 
analyzing the equation satisfied by the error <f> (r, z) := u (r, z) — H ( t ), we find 
that for ro large, around (ro, / (ro)), the C 2 norm of </> is small. 

Observe that u (r, / (r)) = 0. Hence using the definition of <f>, 


H 


( f (r) - [/ (ro) + f (t-q) (r - 

V ^1 + W (ro )) 2 



+ <t>(r,f(r)) = 0. 


(15) 


Differentiate this equation with respect to r, using the fact that d z <f> is small, we 
get | f (r) — f (ro)| is small, which we already know. Differentiate (fl5l) twice, 
we obtain 


H' (•) /" (r) + H" (•) (/' (r) - f (r 0 )) 2 + d 2 r <t> + 2 d r d z <t>f + d z fif = 0- 


Using the fact that d z 4> is small, we deduce /" (r 0 ) is small. 

Indeed, one could continue this process and show that the higher order 
derivatives are also small. ■ 
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Since |/"| is small, the Fermi coordinate actually will be well defined in 
a very large tubular neighborhood of A f u . For later purposes, we need to be 
more precise in describing the size of the region where Fermi coordinate is well 
defined. 

For each point Z = ( r,z ) £ E, we use dist(Z,J\f u ) to denote the distance 
between Z and the curve M u . That is, 

dist (Z,Af u ) = min{| Z — p\ ,p £ M u } • 

Let 7 T (Z) be the set of points which realize this distance. Hence 

?r ( Z ) = {p £ Af u : \Z — p\ = dist (Z, J\f u )} ■ 

For each n, let us consider the set 

D ri := {Z:n(Z)={(r 1 ,f(r 1 ))}}. 

Note that by the triangle inequality, D ri is connected. Then 

D ri = {X (n,t) : t £ (—di (n) ,d 2 (ri))} , 

for some functions d\ and c? 2 (di may not be differentiable). 

Define 

d(n) = min{di (n) ,d 2 in) ,3/ (ri)} 

and 

D ri = {X (n,i) : t £ (-d(ri) ,d(n))} . 

Here the constant 3 has not particular importance. Fix a large constant rg. Let 
B u := U ri>ro D ri . It is worth to be pointed out that in principle the boundary of 
B u could be quite complicated. Let 77 , whose existence is related to assumption 
0 below, be a smooth cutoff function equals 0 outside U r . 1 > T . 0 _iD T . 1 and 

X*r](n,zi) = 1, for Z! £ (-d(n) + l,d(n) - l) ,n > r 0 . 

Similarly, let B+ = H u flE + and we dehne a similar cutoff function rj + supported 
in H^"(note that the notation r) + does not mean the positive part of rf). 

Since B u is not the whole r-z plane, we shall use the cutoff function 77 to 
define a smooth function "Hi by the formula 

Ui = rj H i + (1 — 77 ) 

Here the function Hi is defined through 

X*H x (n, zi) =H(n-h (zi)), 

where h is a small function to be determined later. We emphasize that since we 
are only interested in the behavior of u outside of a bounded set, Hi should be 
regarded as a function only defined on the set E, where 

E = E\|(r,z) : -f(r 0 ) + J7J—j (r - r 0 ) < z < f (r 0 ) - (r-r 0 )j . 
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For any function 0, let 0 s (r, z) = 0 (r, —z ). 

With all the previous notations introduced, we now define an approximate 
solution u as 

u = n 1 + n s 1 + i. 

The function u is defined on the set E, rather than the whole plane E. Let 
(f> = u — u be the difference between the true solution u and the approximate 
solution u. Note that since u is even with respect to the z variable, <f> is also 
even. Introduce the function H[ by 


X*H' 1 (r 1 ,zi) = H'(z 1 -h{r 1 )). 

Then the small function h appeared in the definition of the function Hi is 
required to satisfy the following orthogonality condition: 

[ X* (<) dzi = 0, n > r 0 . (16) 

Jr 

The existence of h is guaranteed by the following: 

Lemma 9 There exists a small function h(in C 2 sense) satisfying (1161) . 

Proof. This follows from similar arguments as that of [19j. The basic idea is to 
use the fact that </> is small and apply the implicit function theorem. We omit 
the details. ■ 

One advantage of using the Fermi coordinate with respect to the nodal curve 
is that h could be estimate in terms of (j> and /. To see this, let N s (u) be the 
graph of the function z = — / (r). For Z = (r, z) € E, let 

V (r, z) = dist (Z,Af (u)) + dist (Z,J\f s (u )), 

that is, the sum of the distance of Z to the nodal line in the upper plane and 
lower plane. Set D (n) = V (n, f (ri)). 

Lemma 10 For r± > r 0 , 

\h (ri)| < C \X*(j> (ri, 0)| + Ce~'^ 2D ^ ri ' > . 

Proof. Since (f) = u — u, around J\f u we have 

X*(j) (n, Zl ) = X*u 0 rx, Zl ) - [H (zx - h (n)) + X*H{ + 1]. (17) 

Setting z\ = 0 in the above equation, using the fact that X*u(r i, 0) = 0, we 
find 

x*(j> (n, 0) + h {-h (n)) +1 + X*U{ (n, 0) = 0. 

On the other hand, by the definition of "Hf and the asymptotic behavior of H, 

ll + X'Hl (rx, 0)| < Ce~'^ 2D< ' ri \ 
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Therefore, 


IMn)| < C|^*<A(n,o)| +Ce~^ D(ri \ 

This completes the proof. ■ 

Following the above arguments and differentiate equation (ED with respect 
to 7 *i and setting z 1 = 0 in the obtained equation, we get 

\ti (n)l < C \d ri x*4> (n, 0)1 + Ce~^ D ^\ 

Similar arguments yield estimates for the higher order derivatives and Holder 
norms. 

One of the main ingredients in the proof of Proposition [4] is to get suitable 
decay estimate for the function <j>. In this respect, we shall first of all prove the 
following estimate. 

Proposition 11 For r± > ro, the function <j> satisfies 
\\X*fi (n, -)|| L =o < Ce ~ r1 + Ce~ D ^ 

r i 

We shall first of all prove this proposition under an additional assumption 
on the size of the Fermi coordinate. More precisely, at this stage, we assume 
that 

i _ 9 n (r) 

\d'\ < -,|d"| < C and d(r) > —(A) 
z o 

Later we will indicate the necessary modification of the proof if this assumption 
is not a priori satisfied. 

Clearly <f> satisfies 

L(j):=-A(/)+ (3w 2 -l) fi = E{u) + P (</>). (18) 

Here the notation E (u) stands for 

u rr + r~ 1 u r + u zz + u — ii 3 , 

which is the error of the approximate solution u, and 

P (fi) = -3 ufi 2 - fi 3 

is a higher order perturbation term. 

For any function 0, we define the projection of 0 onto rj + H[ as 

all_ (07? + .ffi) (fei , , 

f R X* ((v+fPi 2 ) dz 1 

Since the solutions are even in the z variable, we also define 

©II = 0j[ + 
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Finally, set 


0- 1 ^e-e 11 . 

Equation (1181) could be written in the form 

L<j>=[E(u)\ L + [E(u)f +P {<!>). 

As we will see later, [E (u)]^ is small(in suitable norm) compared to </>, and </> is 
then essentially controlled by [E (u)]^ . As a crucial step towards the proof of 
Proposition [TT) we shall take the task of analyzing E (u). By definition 

u = Hi+H\ + l. 

A simple manipulation leads to the following expansion 

E (u) = E (Hi) + E (H{) + 3 (' H{ + 1) (U\ - l) + 3 (Hi + 1) (H{ + l) 2 . (19) 

Certainly the function E (u) is also even in the z variable. One expects that in 
the upper plane E + , the main order of E (u) should be E (Hi) and one of the 
interaction term 3 (Hf + 1) (H\ — l) . 

We first analyze the projection of E (Hi) onto rj + H' 1 . For technical reasons, 
we introduce a small perturbation of / due to the presence of the modulation 
function h. Let p = f + h, where h := y/l + f' 2 h. Throughout the paper, we 
set Co = J R H' 2 (t) dt. The notation O (g) represents a function in n such that 
\0(g)\<C\g(n)\. 

Lemma 12 For rq > r 0 , 


f X* ( V + H[E (HJ) dzi = (l + o (1)) ^ + O ( h"f") + O (ti 2 ) 

JR r l \ a/ 1 + P ) 

+0 (h'f") + O ( h’r- 1 ) + O f e ~ 2V2 ^ r i)) + O (h'h") + O (hf") 

+ O (hf") + O (f" 3 ) + O . 

Proof. In the region where the cutoff function rj + f 0, E (Hi) could be ex¬ 
pressed in terms of the Fermi coordinate (ri, zf : 


E(Hi) = \ r , z) Hi 
= A -^ 2 
d z ,A 


l d r H i + H[ - Hf 


dfHi + dfHi 


+ 


2 A 


dziHl ~ ^ 9riHl 


+ r 1 (d ri Hid r ri+d Zl Hid r zi) + Hi-Hf 


Obviously d Zl Hi = H ', d'f H\ = H". and 

d ri Hi = -h'H',dfHi = -h"H' + h ,2 H", 
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where H' and H" are evaluated at z\ — h (?’i). Since H" + H — H 3 = 0, we 
obtain 

E ('Hi) = A- 1 (- ti'H' + h ,2 H") + (^j- + r~ 1 d r z\ S j H' 

+ 2 A 2 r 1 ^ rr ’ 1 ^ ( _ h'i/'). (20) 


It follows that 

f X* (n+HiEiHi)) dz ! 

*/R 

= [ ri + *A- 1 (-h"H' + h ,2 H")H'+ [ »7 + * + r-1 $-3i 

Jr Jr V 


H‘ 


/ 2 


= -/i" / rj + *A^ 1 H' 2 + / 77 


+* ^zi 


A 


2 A 


+ r i 9 r 2:i I/' 


r/2 


+ 


O (/i' 2 ) + O {h'f") + O (, h'f") + O (h'r -1 ) + o (e~ 2V2 ^j . 


The appearing of the term O 2 '/2 <i ^ is due to the fact that in the upper 

boundary of B u , E (Hi) = O (e"^) . 

To proceed, we shall calculate the second term in the above expression, which 
roughly speaking should be the main order term of the projection. We have 


d Zl A _! 
~2A~ r drZl 


y/rn^ 121 (i +f 2 Y 
A 


r 


( n “ Zi 7 ^f) ' yTT7 


r 


z 1 - 


r 


(1 + f 2 y (! + /■ 


1 + 2z\- 


r 


r 


n V 1 + f 2 

f" 


1 + Z\ 


r 


n^/TTr 

zif" 2 


o 


(1 + f ,2 Y 2 


/2 


o(/" 2 ) 


a +n* a +n 
r zif 2 


3 + 0 (/" 3 ) 
f 3 


nVT+r r 2 (l + n 


o 
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Using this estimate, we get 


X* ( V +H[E(H 1 )) 
c 0 h" 


( 21 ) 


c 0 r 


Co/' 


■ " (1 + ° (1)) 1 .TTF + + riy/l + T* 

+ O (ti 2 ) + O (h'f") + O (.h'f") + O (h'r _1 ) 

+ o (. ti'f") + o ( e - 2V5J < ri >). 

In this expression, we are mainly interested in the first term. We calculate 

rih" , n f" , /' 


+ 


1 + /' 2 (l + ff yfi +1 


72 




1 + /' 2 IVT+7 2 

nf" , y 


2 ti 


1 


1 


, y/l + /'• 


\/l + /' 


(i + /' 2 ) 5 vi + r 


rip 


\/l + i 


+ 0 (/i'r 1 /i") + 0 (/i'r 1 /'') 


+ O (hrif") + O (hrif") + O (h'). 


The conclusion of the lennna then follows from this estimate. ■ 

With the projection being understood, we proceed to analyze the orthogonal 
part. 


Lemma 13 For r\ > ro, the following estimate is true: 

IE mo i! = - (* - *>+ Tprrn) H ' 

+ O (. h"f") + O (e“' /2£> ) + O (h'r- 1 ) + O (h" 2 ) 

+ O (f" 3 ) +o( ?p\ + O (h ,2 ) + O (h'f) + O (h'f ). 

( 22 ) 
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Proof. Let us consider typical terms appeared in (TSUI) . 


A~ x h!'H' - [A~ l h"H'}\ 


= A~ l h"H' - 
h"H' 




1 - 


L(v + *H'Y 

rj + 


■rj + *H' 


r,+*H 12 


0(f"h") 


Note that 


Hence 


i r (v + *h ') 2 

= 0{h")( f (7 f*H') 2 -r, + * [ V +*H'Ah' + 0(f"h"). 
\Jr Jr J 

f f (v+*H') 2 - 77 +* f rj+*H’A H’ = 0 (e _V5/ ) . 

\J M «/ M / 


A~ x h!'H' - [A~ 1 h"H'] I 1 = O ( h " 2 ) + O ( e ~' /2I) ) + O (f"h"). 


Next, we have calculated that 


, -u /" 

_ + r 3 r2l =- 


Zif 


//2 


(l + /'2)2 (l + /«)- 
/' ~H/' 2 


nVl + /' 2 rf^ + f' 2 ) 

+o(r s )+o(^y 


Using this expansion, we find 


+ r ±d r Z l ) H ' 


H' 


%r +r " 15 ^ 1 ) H '~ 

u\ ( f ”' 2 , J 

^ h \{l + f' 2 f + rl(l + f' 2 )j 
+ O(/" 3 )+o(0+o(e-^). 

The rest of the terms always contains small order terms of h! or h and the 
conclusion of the lemma readily follows. ■ 

The next result gives us an estimate for the main interaction term between 
7i\ and Ti\ appeared in E (u ). Let 


Cl = 3 y/2 f H' 2 (s) e' /2s ds. 

Jr 
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Lemma 14 For rq > r 0 , 


[ X* [3 V + H[ (HI + 1) (Hi - 1)] d Zl = - Cl (l + o (1)) e~^ D{ri \ 
JR 

Proof. Since H" = H 3 — iiJ, H\ — 1 = —y/2H[, it follows that 

f X* [3 n + H[ (H[ + 1) («? - 1)] dz! 

Jr 

= -3y/2 f X* [ V + H( 2 (H{ + 1)] dz 1 + o ( e-^^j 

J R 

= -3V2 [ X* (v + H[ 2 ) e-^^dz! 

Jr 

+ o(^J x* (j 1 + H' 2 ) e~ 2V2 ^dzJj . 

Since by definition \z?\ = V (r, z) — \z\ \ , we find 

f X* [3 V + H( (H\ + 1) (H\ - 1)] dz 1 

Jr 

= -3v^ [ X* (rj + H[ 2 ) e" v/5(A '* 2:, - |zi|) (fei 
Jr 

+ O (J X* (v + H[ 2 ) e - 2 %/ 2 ( x-D-\ sl \) dz ^j 

= -3V2 [ X* (V + H[ 2 ) g-'/2(A'*D-|z i |)^ 2 i 

Jr 

+ 0 ( e -+2D(rO^ 

= -3V2e~' /2D{ri '> [ X* (H[ 2 ) e^^dzi 

Jr 

+ o(e-^ D W). 

In the last equality, we have used the fact that in E + , 
\X*V(r 1 ,z 1 )-D(r 1 )\=o(z 1 ). 


Remark 15 Clearly the function D (n) is strictly less than 2/ (n). However, 
due to the fact that f is small, 

2f (n) - D (n) = o(f (n)). 

It is exactly this fact that makes it possible for us to analyze the equation satisfied 

by /• 


18 



The previous result analyze directly the error E ( u ) from the definition of u. 
In the next result, we will express the projection of E ( u ) onto r] + H[ in terms 
of the function <j>. The main equation we will use is equation m satisfied by 0- 

Lemma 16 For each n > ?'o, 


X* ( ij + H[E (u)) dz x 

R 

o (llO* (n,OIlL) + o (||0* (n, OIL e-^)) + o (||0* (r 1} OIL II^O* (n, OIL) 

+0 (||3 ri 0* (rr, OIL e-^ /Cri) ) + O (||O ri 0* (r 1; OIlL) + O (110* L, OIL \K<t>* (n, 0||J 
+o (/" \\d zi (/)* (n, OIL) + o {f" ILL* (n, OIL) + o (/"' ILL* L, OIL) 

+o (r- 1 ILL* (n, OIL) + O (r- 1 / 7 ILL* (n, OIL) • 

Proof. Multiplying both sides of (fTSl) with ?/ + U( and integrating in R, we get 

[ X* (r/ + H , 1 E(u))dz 1 = [ X* (ri + H' 1 L(f>) dz\ — [ X* (r] + H[P (0)) d Zl . 

Jr Jr Jr 

Let us calculate the term J R X* (r] + H' 1 L(j)) dz\. By formula (fl3l) of Laplacian 
in the Fermi coordinate, 


/ X* (r] + H[L^) dzi 
Jr 

= [ X*(r,+H[) 


4-1.92 | l dziA d 1 dr i A Q 
^ + 2 A dzi 2 A 2 dri 


<j>* dz\ 


- f X*(r] + H[)r 1 [d ri (l)*drri + d Zl (j)*drZi]dzi 
Jr 

+ j [-X* ( V +H[) d Zi (j>* + X* ( ; n +H[ (3 u 2 - 1) 0)] dz\. 

Jr 

We first estimate f R X* (rj + H[) A~ 1 d 2 1 (f>*dzi. Differentiating the identity 

[ X* (r/+L'0) dz 1 = 0 (23) 

Or 

with respect to ri, we obtain 

[ X* (: r,+H[ ) d ri cj ) *dz 1 = - f d ri [X* ( ri+H[ )] 0*d*r 
Jr Jr 

= - [ d ri (r?+*) H'pdzi + [ rf r * H" h' <f)* dz\. 

Jr Jr 


Here El' and H" are evaluated at z\ — h (ri). Therefore we could apply Lemma 
[HU and obtain 


/ x* (L^RL* 

Jr 

= o (no* (n, OIL + O (110* (rr, OIL ILL* (n, OIL) ■ 
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Similarly, differentiating C3l) with respect to r± twice, and using the fact that 
^ _1 - TTF 1 = ° (/") > we S et 

f X* {i 1 + H[)A~ 1 d 2 r ^*dz 1 
Jr 

= rrr L x * ^ +H[) d2 ^ dZl+0 (ri ’- } lloo \f" wo 
= -TTn L dri [x * (??+H[)] drA * ~ TTF L ^ {n+H ' l)] r 

+ °{\\ d r x 4 >* ( r l, ") Hoc \f" ( r l)l) 

= o (||3 ri 0* (n, OIL e-^ /(ri) ) + o (lIM* (n, OIlL) 

+ o (||0* (n, OIL e-^ /(ri) ) + O (110* (n, OIL 11^0* (n, OIL) 

+ 0(110*^,OIL K0*(r l5 OIL) - 

Next, using d Zl A = O (/"), we infer 

/ %^0*X* (77+ifQ d*i = O (/" ||9 Z1 0* (n, OIL) • 

Jr ^ 

Observe that d ri A = O (|/"| + |/"'|), thus 

jf %^s ri 0*x* (77+^0 d Zl = o (in im0* (n, oil)+o (in 11^0* n, oil) • 

We also have 

[ r _1 X* (v +H i) [<9 ri (f>*d r ri + d Zl (j)*d r Zi\ 

Jr 

= o (r- 1 ||<9 ri 0* (n, OIL) + 0 (r -1 l/'l 11^0* (n, OIL) ■ 

Another term we need to estimate is 

/ x* (nn) C0* + [ x * (v + h[ (3 n 2 - 1) 0) 

Jr Jr 

= f X* ( V +H[) d? x 0* + / ** ( V +H[ (3 H 2 - 1) d>) 

Jr Jr 

+ 3 [ x* (ri+H[ (v 2 - H 2 ) 0) . 

Jr 

To handle them, we integrate by parts for the first term, and for the last term 
we use the fact 

77+ (u 2 - H 2 ) H [=0 (77+ (HI + 1) H[) = O (e~ D ). 

We conclude that 

[ X* (77 +H[) d 2 Zi r + [ X* (77 +H[ (371 2 - 1) 0) 

Jr Jr 

= o (||0* (n, OIL e_v5/ ) + o (110* (n,OIL 11^0* (n, OIL) ■ 
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Finally, since P (</>) is higher order term of (f>, the estimate of the term 
f R X* (r] + H' 1 P (</>)) dzi is trivial. Combining all these estimates, we get the 
desired result. ■ 

With the error term E ( u) being understood, we would like to recall some 
properties of the operator L. The next result are essentially proved in [10], 
although here we need to do slight modifications due to the presence of boundary 
terms. For each r± > ro, let 

K rt ■= e E,-f(r i) + j ^-y (r-n) <z< f(n) - ( r ~ ©)} > 

and J ri := E\K ri . Note that J ro is simply the set E introduced before. 
Lemma 17 Suppose p is a solution of the equation 

Lp = 0 in J ri . 


and p (r, z) = p (r, —z) , p\ = 0. Then 

\p* (s, z)\ < C||^|| ioo ( 9Jr je r ' 1 " s + ||0|| ioo( ' Jr j , for each s > n. 

Proof. Let 9 be a cutoff function such that 

„ , f 1, if Z £ J ri and dist (Z, dJ ri ) > 1, 
d ^ = { 0 ,ZiJ ri . 

Consider the function p := 9p. Then p satisfies 

Lp = 9Q + A 9p + 2 \79Vp. (24) 

Modify the function u in E\J ri if necessary, we could get a function, still denoted 
by u, whose nodal curve are almost horizontal and around its nodal curve, 
it looks like the heteroclinic function H. With slight abuse of notation, the 
corresponding linearized operator will be still write as L. We could also assume 
that p satisfies the orthogonality condition. 

Denote 0i = 00, 02 = A 9p + 2 S79\/p. For * = 1,2, consider the equation 

Lpi = Qi + ki(ri)r] + H[ + (kir] + H[) s (25) 

where pi, ki are unknown functions and we require f R X* {ii + H[pi) dz\ = 0, i = 
1,2. By the results of [10 , one could find solution (pi. kA to this equation. This 
pair of solution is indeed unique and hence we have the following estimate for 

p x : 

H^ilLoo < C , ||0 1 || L » . 

For the function p 2 , since 02 has better decay property than 0i, we could 
estimate 

\T 2 (r,z)\<C\\e 2 \\ L „ e r '~ r . 
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Notice that 

L (Cpi + (p 2 ) = @1 + ©2 + (ki + fe) H'i + [(fci + fe) r] + H (] . 

But the equation 

LCpj, = 0! + 0 2 + k 3 r) + H[ + [fc 3 77 + i?(] s 

also has a unique pair of solution (<p 3 ,k 3 ) if we require f R X* (i^ + H[(p 3 ) = 0. 
This implies that (p = Cp\ + </? 2 , the conclusion of this proposition readily follows. 

■ 

Once we have an estimate in L°° norm, we could get a corresponding Holder 
norm estimate by Schauder’s elliptic estimate. Now we proceed to the proof of 
Proposition [TT] 

Proof of Proposition illl Recall that 

/ X* {r, + H[E(u)) 

Jr 

= [ X* (n +H[E (Hi)) + f X* (r, + H [3 {U{ + 1) (%\ - l)) 

Jr Jr 


+ 


[ X* ( V +H[E (%{))+ [ X * (r,+H(3(Hi + l)(Hf + l) 2 ) . 
Jr Jr v ' 


We have already analyzed the first two terms. For the third term, slightly abuse 
the notation, we have 

[ X* (j] + H[E (HD) dz\ = [ X* ( 1 n +H [) H -1 (-ti'H' + h' 2 H") dz x 
Jr Jr 

+ [ X* ( V +H[) 

Jr 


+ r 1 d r zJ S j H'dzi 

+ J^X* (r fH[) - r- 1 ^ (tiH') dzt. 

It should be pointed out that here A, h are evaluated at r 2 and H' is evaluated 
at Z2 — h (r 2 ) , and ti means the derivative with respect to 7-2((r2, Z2) being the 
Fermi coordinate with respect to z = —f (r)), rather than r\. This indeed makes 
the analysis a little more complicated. We have 

[ X* ( ri + H [) A~ x (-ti’H’ + ti 2 H") dz x 
Jr 

+ f X* (r,+H[) r-'drr 2 ) (tiH') dz 1 

Jr 


2 A 2 


= o 

For the same reason 


(e-^ D ) . 


J X*(j] + H() 1 9rZ 2 ^j H'dzi = o(e D ). 
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It follows that 


/ A'* ( r, + H[E (HI)) dz x = o (e~' nD '\ . 
Jr ' ' 


Similarly, 


J X* (v + H[3 (Hi + 1) (H{ + l) 2 ) = o (e-^) . 
It then follows from Lemma [T2] and Lemma [TJ] that 

[ x* (n + h[e(u)) 


= -(l + o(l)) 


c 0 h" 


+ 


Co/" 


+ 


Co/' 


_1 + /' 2 (l + /'2)2 r ly /i + ff 2 
- Cl (1 + o (1)) e _v ^ D ( ri ) + O ( ti 2 ) + o (/z'/") 

+ O (tif") + O (h'r~ l ) + O ( ti'f"). 

This together with Lemma [TG] leads to the following equation: 

c 0 /' 


- (1 + o(l))cie 


-y/2 D 


(26) 


cp f" Cp h" 

(1 + f ,2 f 1 + /' 2 riy/l + f 2 

= O (ti'f") + O (ti 2 ) + o (tif") + O (tir~ l ) + o (tif") 

+ o (ll/>* (n, OIlL) + o (|| 0 * (n, OIL e-^)) + O (||</>* (n, OIL \\d ri f (n, OIL) 

+ o(ll^* (n,OILe- V5/(ri) )+o(l|a ri 0* (n,OllL)+o(ll^ fa, OIL K** fa,OIL) 
+ o (/" um* (n, OIL) + o (/" ||s ri «/>* (n, OIL) + 0 (/'" W 9 ^* fa, OIL) 

+ O (r" 1 ||9 ri ^* (n, OIL) + ° ( r_1 /' fa, OIL) • 

On the other hand, 

Lf=[E(u)} ± + [E(u)] ll +P(f). 

By Lemma [TCI [E (u)]^ + P (f) is small compared to <j>. In view of Lemma [T71 
we need to estimate [E (m)L . 

By m , 

[E (ti)]- 1 - = (E (Hi)} 1 ' + [3 (H{ + 1) (H\ - 1)] X 
+ [E (Ht)} 1 - + [3 (Hi + 1) (H{ + l) 2 


Let us analyze the term [E (H{)\ . We know that 


E (H{) = A -1 (-ti'H' + tiH") + 


d Z2 A 
2 A 


+ r 1 d r z 2 ) H' 


d r2 A 
2 A 2 


— r l d r r 2 ) tiH'. 
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Hence in the upper plane, due to the presence of H' in each term, we find that 
directly that 

[EiHl)} 1 - = o(e~ V2f ^ . (27) 

Note that this is not an optimal estimate and each term in E (Hf) is mulitplied 
by h ', h", f" or r^ 1 /'. However, these terms are in general not evaluated at r±. 
One also has 

3 (Hi + 1) (Hf + l) 2 

The above estimates combined with Lemma QjJ] and [Uj] leads to 



[E (u)] 1 - = [E (Hf )] X + O {f" 2 ) + O (r" 2 /' 2 ) + O (h"f") + O (/iV" 1 ) (28) 

+ O ( h " 2 ) + O ( h ,2 ) + O ( h'f" ) + O ( h'f ’") + O . 

Notice that there is a term f " 2 appeared in the right hand side. Insert (17B1) 
into (l2Sl) we get 


(Eiu)) 1 - 


= [E (Hf)] X + O ( r~ 2 f 2 ) + O (h'r- 1 ) 

+ O ( h ,2 ) + O {h'f"') + O 
+ O {h'h") + O (. hf") + O ( h" 2 ) 


(ll<f>* (ri,-)llL)+0( 

II d ri r (n.-)liy 

(ll d Zl <t>* (n,-)llL)- 



(29) 


Note that the term involving h and its derivatives are essentially controlled by 
ct>. Using similar arguments as that of Lemma 1171 we find that for s > r \, 


\<t>*(s, Zl )\<C\\e D \\ L ~ i{ri , +oo)) +C\\r 2 f' 2 \\ Laoari+oo)) + M L ~(dj ri ) eri '■ 

(30) 

We remark that there is the term [E ('Hf)] appearing in (17iH) . Initially, by (1271) . 
it is only of the order O (e — ■') . But we would using a bootstrap argument to 
show that this term is indeed of the order O (e~ D ) ■ 

The estimate could be applied repeatedly. This yields that for some 
large but fixed constant l, 


C 


4 (r, z) < Ce - D (’■-0 + —, + C ||</>|| 


< Ce~ D ^- l) + 


C 

{f-lf 


\L™(dJ r -l) ' 


c 2 


D{r 2l) + 77-7n2 + ll^llL“(9J r _ 21 ) e “ 


(r - 21) 


-l 


c 
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The last inequality follows from 


k 


E 


(r - jl ) 


o~i l 



and by Lemma[Gl /' is small, which together with D = (2 + o (1)) / implies that 


Y^ e - D (r-jl) e -jl < Ce~ D{r \ 

j =i 

This concludes the proof. ■ 

Up to now, we have assumed 0) . But a priori, we don’t know whether this 
is fulfilled or not. It the sequel, we sketch the proof without this assumption. 
The basic principle is, using equation (l26l) and similar arguments above, one 
could prove that assumption CT) indeed holds. 

Proof of Proposition !!!! without Assumption (A). The main reason that 
we introduce the Assumption 0 is that in the the error E (u) is related to 
e^ 2v/2c * and we could control it by e~ D . Now without Assumption 0. to handle 
the term O , we would like to use the following basic geometrical fact 

concerning the Fermi coordinate (see m for a proof): Using the fact that /' (rq) 
is small, there exists a point f(f may be equal to ri) with 

\r-ri\ = o(min{di (rq) ,d 2 (rq)}) 


such that 


min{cii (rq) ,d 2 (rq)} > 


1 


(31) 


I /"(r)r 

Let us still assume |cZ'| < \ and |J"| < C for the moment. Then repeating 
the previous arguments, one could show that for r\ > rg, 


If" (rq)| < Ce~ ri + — + Ce^ 0 ^ + C 

r\ 


-2V2 d 


L°°(r i,+oo) 


We claim there exists a universal constant Co such that 

d(rq) = 3/ (rq) for rq > C 0 . 


(32) 


(33) 


Suppose to the contrary that m is not satisfied at some large point s. 


Without loss of generality we assume 
we could find Si with 


,-2\/2 d 


0 -2V2d(s) 


L°°((s,+oo)) 


. By dSH) 


|si - s| = o (min {di (s), d 2 (s)}) = o (/ (s)) 

such that 

min {di ( s),d 2 (s)} > „ . 


(34) 
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We would like to show that for this point si, necessarily 


1 

LfM 


> 3/(s). 


( 35 ) 


Once this is proved, we then get a contradiction and hence the proof will be 
complished. 

To prove (Idl^l , we could assume Si ^ s, otherwise using estimate (l32l) . we 
obtain (l35l) . For the point s 1 , we have 

I/" (si)| < Ce~ Sl + — + Ce~ V * DM + C 

si 

= Ce~ Sl + — + Ce~ V2D{si) + Ce- 2V ^ S{Sl) 

si 


^—2y/2d 


Tj°° (fit 1 .-l-o 


for some point si > si. If Ce~ 2V2d ^ < \ \f" (si)|, then 
I/" (ai)| < Ce~ Sl + — + Ce~' /2D ( S1 \ 

Sl 

which implies (131^1 . Hence we could assume 

Ce-^dis,) > I|/"( Sl )|. (36) 

This in particular implies that 

For the point si, by m , we could find a point S 2 with 


S 2 > si + o (min {di (si), d 2 (si)}) > Si - In / ( s ), 


such that 1 

min {di (si),d 2 (si)} > 

Note that by (13(1 , 

Ce-7^T >i|/"( Sl )|. 

Now repeat this argument and we get a sequence {s ra } with 


^n+l ^ s n In 


1 

ww 


(37) 


such that 
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If the process stops at the step no < / (s). Then 

/" (s no ) < Ce~ Sn ° + — + e~ D ( s "o). 

Sn 0 

By (1571) , s no > s — f (s) and D (s„ 0 ) > / (s) . This then implies (1551) . On the 
other hand, if the process doesn’t stop at the step no > / (s), then the sequence 

pn ~ w &\ satisfies 

p n > Ce Pn+1 . ( 38 ) 

(1551) leads to the inequality 


Pi > Cnl > Cf 2 (s), 


This also implies (1551) . 

Let us we return back to the assumption that | d' | < ^ and | d"\ < C. Essen¬ 
tially, this assumption is to make sure that the cutoff function rj has bounded 
first and second derivatives and thus the error of approximate solution could be 
uniformly controlled. However, with the original definition of d this assumption 
may not be true. To overcome this difficulty, we should modify the function 
d(that is, one needs to modify the domain B u ). Let us be more precise. Fix a 
small positive constant S. We define new function d by 

d (ri) := inf {d(s) + 6 (ri — s) : r 0 < s < n; d(s) — 5 (r\ — s) : s > n} . 

Modify d such that it becomes C 2 . We then define the domain B u using this 
new function d. Using similar arguments as before, we could show (1331) . 

Now m implies that the size of the Fermi coordinate is actually large 
enough for our purpose, that is, Assumption © is indeed satisfied. ■ 

With all these preparation, we proceed to prove Proposition [Tj the main 
result of this section. 

Proof of Proposition [4] and Theorem [2j We split the proof into several 
steps. 

Step 1. We first show that there exists a universal constant C, <5 > 0 such 
that 

p{r) > C + + <5^ In r. 

We would like to use equation (1551) . Note that there is a term involving the 
third derivative of / in this equation, although one expects that f" decays like 
r^ 3 , a priori we don’t have any decay information for it. However, we at least 
know that f" (t\) tends to zero as r\ tending to infinity. Our aim is to show 
the following estimate(not optimal): 

I/'" (r-i)l < CrC 1 + Ce~ D ^\ (39) 

To obtain this estimate, we would like to differentiate equation (l26l) with respect 
to n. Using the fact that and /^ are small for n large(at this stage, we 
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don’t know the decay rate for these two terms, but they are multiplied by terms 
related to <j >), we obtain the following estimate (which is not optimal but enough 
for our purpose) 


/"' + v / T+7 *h'" 

= O (h") + O (rq 1 ) + o 

+ O ( h ') + O (\\d>* (n, OIL) + O (ll d ri f (n, OIL 
+ O (\\d 2 ri r (n, Oil J + O (WdrM* (n, OIL) ■ 

Hence by the C 2,a estimate of 0) to prove (l39l) , it will be suffice to obtain an 
estimate for the function h!". which is essentially controlled by df . 1 (j>* . To achieve 
this, we consider the equation 

L<j>=[E(u)] ± + [E(u)] 11 +P(</>). ( 40 ) 

Here we have in mind that by Lemma [Tbl [E (u)]^ is expressed as a small order 
term of <f>. We also have the estimate (l29l) for the source term [E (u)] 1 ". Now 
differentiate equation (HOI) with respect to rq, use the L°° norm estimate of </> 
and the fact that 


d ri [E (u)] 1 - = O (L 1 ) + O ( e ~ D ) + o (h "'), 


we find that for rq large, 

\d^(j)*\ < Cr^ 1 +Ce~ D . 

This then implies d39l) . 

With the decay estimate of /'" available, equation (l26l) could be refined to 


( np' \ 


—ri (l + o (1)) e^ D + O (rr 2 ) := h + h- 
Co 


(41) 


Here p = f + \J\ + f' 2 h is the function introduced in Lemma [12] At this stage, 
one of the technical difficulty is that this equation involves the function D. As 
we mentioned before, we expect that D is very close to 2 p. But without a priori 
estimate for /' and f,\D — 2p\ in principle could be large. In any case, we at 
least know that D < 2 p. 

To proceed, we shall fix a small constant S > 0 which will be determined 
later. For any interval (ti, t?) C (?*o, +oo), integrating equation (UTl) in n yields 

hp' (t 2 ) _ hp' (H) 

V 1 +P’ 2 (*2) L 1 +P' 2 ( f i) 

Keep in mind that J)* 2 1\ is always positive, and by \I% (r)| < Cqv~ 2 , 




< 


Co 

tl 


(43) 
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Set t* := max {^f 2 -, ro} . We consider two cases. Case 1: 


t*p' (f) ^ V2 

VI +P ,2 (t*) ~ 2 

Then by (1T21) and (H51) . for all t > t*, 

tp' (t) > V2 

V 1 + p ,2 (t) ~ 2 ' 

Case 2: _ 

t*p' (f) V2 

a/ 1 + p' 2 (t*) 2 

In this case, let be an interval such that 


6 . 


6 

2 

<5. 


V ft) 

V 1 + P /2 (i) 



+ 2<5,Vf £ (f*,f) . 


Then the fact that \p'\ = o(l) and a simple integration yield 


P(t) < 



In t + C, Vi £ (f*, f) . 


This in particular implies that D ( r ) — 2/ (r) = o (1) and hence for any ti,t 2 £ 

C 


rt2 se-^ds> ° 


.& V 28 , 6\/28 ' 

l l Z 2 

We choose 5 such that 6\/2 5 < 1. Then by inequality (T43l) , if t 2 is large, 


pt 2 


h + 


h > o. 


This implies 

t2p' (h) ^ bp' (ti) 

V 1 +p' 2 fe) V 1 + p' 2 (*i) 

Combing the analysis for Case 1 and Case 2, we find that there exist universal 
constants <5 and C\ such that 


/2 

rp' (r) > — + <5, r > Ci. 


( 44 ) 


This proves Theorem [5] 

With the lower bound (Hil) available, one could show that 

e~^ D = e~ 2V * p + o(r-V +a A . 
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It also follows that 


p' (r) = —I - O (r 
r 


p'(r) = -+0 (r- 1 -*) , 


k 


for some a > 0. Here k is the growth rate of u. Therefore, 


p (r) = k In r + O (r “) + C. 


With this estimate at hand, we find that u n converges to a two-end solution uq 
strongly. ■ 

3 Moduli space theory of two-end solutions 

3.1 Preliminary results 

Generally speaking, the structure of the set of bounded entire solutions to the 
Allen-Calm equation could potentially be very complicated. However, if one 
impose certain conditions at infinity for the solution, then it could be simpler. 
The moduli space theory for multiple-end solutions of the Allen-Cahn equation 
in R 2 was developed in [9j. This theory tells us that if u is a 2fc-end solutions 
in R 2 and nondegenerate, then around u (in suitable sense), the set of 2fc-end 
solutions is actually a 2fc-dimensional manifold. This fact has been used in an 
essential way in the classification of four-end solutions of Allen-Cahn equation 
in R 2 ([193). In this section, we would like to develop the corresponding moduli 
space theory for two-end solutions in R 3 . Our main result states that the moduli 
space of two-end solutions has a structure of real analytic variety of formal 
dimension 1. 

To begin with, let us recall some preliminary results about real analytic oper¬ 
ators. Compared to C°° operators, real analytic operators has better structures. 
We refer to 0, 0, 0 for more details. 

Let X and Y be Banach spaces and U an open subset of X. We first recall 
the notion of real analytic operator. 

Definition 18 A map F : U —> Y is real analytic at Xq € U if there exists a 
S > 0 such that 



where nik is a symmetric k— linear operator, and there exists r > 0 such that 


SUp r k ||77T,fc || < + 00 . 
fc>o 


The function is said to be real analytic on U if it is real analytic at every point 


ofU. 
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Let F : U — > Y be a real analytic functional. Suppose that dF (x ) is a 
Fredholm operator of index 1. Assume there exists a map A : (0,e) — > Y such 
that F (A (s)) = 0, and dF (A (s)) : X — > Y , is surjective for all s € (0, e). Let 

S = {x G U : F (x) = 0}. 

The following theorem has been proved in the book of Buffoni and Tolarid [5]. 

Theorem 19 Suppose all bounded closed subsets of S are compact in X. Then 
there exists an extension of A, denoted by A : 

(0, +oo) —■> X , 

satisfying: (1) A is continuous. (2) The set of points where dF is not surjective 
has no accumulation points. (3) One of the following happens: (i) ||A(s)|| —» 
+oo, as s —► +oo; (ii) A (s) approaches the boundary ofU as s tends to infinity; 
(Hi) A((0,+oo)) is a closed loop. 

Basically, this theorem tells us that if one has a real analytic variety which 
comes from the zero set of a real analytic operator and its formal dimension is 
one, and assume further that on the variety there are some points where the 
operator is surjective, then under certain compactness assumption, one could 
find a continuous path of solutions where at most finitely many solutions are 
degenerate (the linearized operator is not surjective). 

3.2 The real analytic structure of the moduli space 

Let u be a two-end solution. Then u satisfies © and there are constants k and 
Cfc € R such that 

IkOv) — H(- - fchrr — c fc )|| LOO([0 +oo)) -> 0, as r -> +oo. (45) 

We also assume k > y/2. The moduli space theory of noncompact geometric 
objects with controlled geometry at infinity (minimal surfaces with finite total 
curvature, singular Yamabe metrics, constant mean curvature surfaces with De¬ 
launay ends) has been developed in [2Tj, [23]) [25]. In this paper, we will not 
investigate the general moduli space theory for the Allcn-Cahn equation in di¬ 
mension three. Instead, we shall only consider those solutions satisfying ©■ 
Our aim is to show that the set of solutions of © satisfying (14511 has the struc¬ 
ture of a real analytic variety with formal dimension 1. Additionally, if a solution 
u is nondegenerate, then around it, this real analytic variety is actually a one 
dimensional real analytic manifold. 

Consider the function 

fk,b (0 = A; cosh -1 (fc -1 r) + b. 

where k and b are real parameters, k > \[2. Obviously, / represents a verti¬ 
cally translated catenoidal end and we have in mind that / is the asymptotic 
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curve of the nodal line of a solution u. The moduli space theory for two-end 
solutions will, roughly speaking, state that around a given solution u there is a 
one dimensional family of solutions, with k or b being possible candidates for 
the local parameters. To make this statement more precise, we need to analyze 
the mapping property of the linearized Allen-Cahn operator around u. 

Let us introduce some notations. As in Section [2j we shall also carry out 
the analysis through the Fermi coordinate. Let B u be the domain where the 
Fermi coordinate of z = f (r) is well defined(suitably modified if necessary). We 
still use 77 to denote a smooth cutoff function supported in B u . Similarly, we 
have the cutoff function rj + supported in B u D E + . Similarly, we have the map 
X : (j"i, z\) -> (r, z). 

We need to work in suitable functional spaces Si , S 2 which will be described 
now. Let C\ > 0 be a fixed constant and 5 > 0 small. A C 2,a function S £ Si, 
if and only if S (r, z) = S (r, —z) and it satisfies the following condition: In E + , 
when n > Ci, 


X *S (r lt zi) = p (ri) X*r] + (n, zi) H' ( 21 ) + X*ip (n, 21 ), 
where f R X* {p + il’) H'dzi = 0, 


(1 + n) p 


and 


(l + n) 3 p' 


(1 + nfx*^ 


C 2 .' 


(1 + Tlf p" 


< + 00 . 


c °.■ 


< + 00 , 


Similarly, a C 0,a function S £ S 2 if and only if H (r, z) = S (r, —z) and it 
satisfies the following condition: In E + , when r\ > Ci, 


x*e (n, zi) = p (n) x*r] + (n, z x ) H' (zi) + x*ip (n, zi ), 
where f R X* (p + il’) H'dzi = 0, 


(1 + ri) p 


C°>' 


< + 00 . 


and 


(1 + r i ) 4 X*ip 


C°-‘ 


< + 00 . 


Let £ be the linearized operator of the Allen-Cahn equation around a two- 
end solution u., that is, 


£ = C u := + r 1 d r + 1 — 3 u 2 . 

For each point p, let (ri (p ), z\ (p)) be its Fermi coordinate with respect to 
fk,b- Let p' s t be the point whose Fermi coordinate with respect to the curve 
2 = fk+ s ,b+t is still equal to (n (p) , zi ( p )). For p in the upper half plane, now 
let $ Sl t be the map defined by 

®s,t ip) = V + p's,t + (! - ^ + ) 
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while for p = (r, z) in the lower half plane, we let 


(p) = -$ s ,t ( r,-z ). 

In this way, we have defined the family of maps <f> st which is even with respect to 
the r axis. For |s|, \t\ small, is a diffeomorphism. By definition <f>o,o (p) = P- 
Similarly, we introduce the family of diffeomorphism 

®a,t (p) ■= PP^,t + (! - V)P- 

Similar arguments as that of Section [I] shows that at far away an axially 
symmetric two-end solution u could be written as: 

u = Hi + H{ + 1 + <t>. 


Here 

Hi = pH 1 + (1 - 77) 

I-nil 

where 

X*Hi (n,zx) =H (Zi - h (n)) 

and 

[ X* (r) + (/)) H'dzi = 0. 

J R 

Notice that 4> decays like r~ 4 at infinity. 

For £ R 2 ® Si, we then define a family of function u a ,t,ip such that 

at far away 

u s ,t^ = Uxo 'F-j 1 + [Hi o®’}] s + (<j> + ip) o $-1, 

and in a fixed large ball, u s , t ,^ = u. Clearly, for u 0 , 0,0 = u • Consider the non¬ 
linear map N : R 2 © Si —> S 2 given by 

(S, t, Ip) ->• [A Us,trf + U s , t - u 3 s , t ,^\ o $ s , t . 

The reason that N maps R 2 © Si into S 2 lies in the asymptotic behavior of u, 
see Section [2] We have 

Ss X ,£,■!/> ,t/j] ^ © S s • X [A^ ; 

dtX (s, t , 'iji) — ^t^s,t,ip\ O *&s,t © ' X [A U s ,t,il; © U s ,t,ij> ^ •> 

d^,N(s,t,ip)G = [Cu sM G\ o<f> M . 

In particular, although DN is not exactly equal to £, it is a small perturbation 
of it for s, t small. Observe that when (s, t, if>) = (0, 0, 0), DN actually is equal 
to C. Let 


71 9 S U S ,t,il;\ (s,t,£/>) = (0,0,0) ) 

72 = dtu s ,t,ip\( s ,t,ip)=(o,o,o)- 
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Introduce the deficiency space 


T> =span { 71 , 72 } • 

The next result concerns the mapping property of the linearized operator C and 
is the main result of this section. 

Proposition 20 The operator £ : Si®D S 2 is a Fredholm operator of index 

1 . 


Proof. Let us prove the result under the additional assumption that / (0) is 
large and H/'H^ , ||/"|| ioo are very small. In the general case, we could modify 
the function u inside a compact set into a function whose nodal lines are almost 
parallel, and use the fact that the corresponding linearized operator is a compact 
perturbation(thus the Fredholm index is preserved) of C. We shall split the proof 
into several steps. 

Step 1. For each function 0 £ S 2 , we shall find a solution S to the equation 


CE = 0. (46) 

By the definition of 52 , in E + , for rq large, 

X*Q (rq, zi) = w (rq) X*r] + H' (zi) + X*tp (rq, 24 ), 

for some function oj and <p, where J R X* (r/ + ip) H' = 0. To solve (l46l) , adopt the 
same notation as in Section [2j it will be suffice to solve the following system 


(47) 


I (£3)f = 0“, 

1 (/S) ± = 0 - L . 

For convenience, we recall the expression of £in the Fermi coordinate. 
r — A~ 1 g 2 -l q 2 I ^ dzi A _ 1 ^1 ^ g 

l, — a a ri + a Zl + 2 a z x 2 ^42 ° Ti 

+ r _1 ( d r r 1 d ri + d r zid Zl ) + 1 - 3 u 2 . 

Let us first compute the action of C on functions of the form £ (jq) rj + H' ( 21 ). 
Using m we get 


c (fv + H') = a~ 1 Cv + h' + 


1 d Tl A 


2 A 2 r (1 + f' 2 ) B 


CV + H' 


1 d zl A 


r 


£r? +H" 


2 A T \J 1 + f' 2 

+ [H'" + (1 - 3 u 2 ) H'] £ 77 + 

+ 2A~ 1 £d ri ri + H' + A- l fdlj t + H' + fd Zl r,+H" + £^, ? + H' 

fd ri r] + H' 


1 d ri A 


1 


+ 


2 A 2 r (1 + f 2 ) B 

r 


1 d Zl A 

2 A 


-yTTT 7 


fd ZlV + H'. 
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The coefficient before £' has the estimate 


1 d ri A 1 

2 A 2 + r (1 + f' 2 ) B 




Next we shall estimate the coefficient before £. Recall that 


1 d Zl A _ f 

2 A rsjl + f' 2 


O (r-) . 


We also have 


H’" + (1 - 3it 2 ) H' = 3 (H 2 - u 2 ) if' 

= 3(if + u)(i?-u) fT. 


Since 


u — H = H (z\ — h (n)) - if (21) + O (r" 2 -“) 

= -h (n) £?' + O (r _2_Q ) = o (r _2_Q ) , 

hence H'" + (l — 3it 2 ) H’ — O (r -2_a ) . As for those terms involving derivatives 
of ?? + , they could be estimated by O (r -2-a ) . Combining all the above estimate, 
we get 

C (&+H 1 ) = A~ 1 C'v + H' + (J- T O (rf 3 )^) Crj+H' + O (r~ 2 ~ a ) f. 
Introduce the operator 


Ki : £f C (£,rj + H r ) ■ n + H'dz 1 . 
Jr 


It will be important to analyze the mapping property of K\. Let us consider the 
equation 

Ki£ = oj. (48) 

This is a second order ODE. Suppose are two linearly independent solu¬ 

tions of the corresponding homogeneous equation: K\ dj = 0. We could assume 
is bounded near 0. Let y be a cutoff function equals 0 in (0,1) and equals to 1 
in (2, Too). Then using the variation of parameter formula, one could show that 
K i is an isomorphism from the space S\ ® span {x'J'i} to <S> 2 . Here a function 
p G S 1 iff p G C 2 ’ q (R) and 


(1 T nfp 


T 


(l+i'i) 3 ?' 


T max 


(1 T n) p r 




< Too, 


and a function p £ <S 2 iff p € C 0,a (R) and 


max 

S 


(1 Tri) 4 p 




< Too. 
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Define the operator 


K 2 :^ [C^+H'+l^+H 1 ] 3 )^ . 

To solve gZD , it suffices to find solution (£, ip) to the system: 

r K 1 £ = w- hv+H'CApdzx, / 4q x 

\ (Cipt =<p~ K a t [ ’ 

Let ipo be the solution of 

{£ipo) ± = ~K 2 (X& 2 ) , 

satisfying f R X* (ip 0 r] + ) H'dzx = 0. Note that at infinity, dx behaves like O (lnr) 
or O (1). Hence K 2 dx behaves like O (r~ 2 ~ a ) . This implies that ipo = O (r -2-ct ) . 
Thus f R X* (Cipor] + ) H'dzi = O (r -4 ) . Thanks to these estimates, one then 
could use a fixed point argument to get a pair of solution (£, ip) for (H51) with 
the form 

f f = CiX'd'2 + { 0 , 

\ ip = C\1po + ip , 

for some constant ci, where £0 £ <Si and ip £ S':.(Note that ip 0 does not belong 
to Si). Then we get a corresponding solution S to (H71) . We emphasize that at 
this stage, we still don’t know whether the solution E belongs to Si, due to the 
non-decaying part 

S := x (ri) d 2 (n) ri + H' (zx) + [x (rx) d 2 (rx) r] + H' (zx)] 5 + tp 0 . (50) 

Step 2. Investigate the homogeneous equation 

CC, = 0. (51) 

Firstly, we wish to find a nontrivial solution to (l5l1) . To achieve this, we use 
the fact that i9i is in the kernel of K \. A perturbation argument could be applied 
similarly as in step 1 to get a function 5 0 solving £5 0 = 0 ,where S 0 is around 
? 7 + $i (ri) H' (zx) + [v + dx (ri) H' (zi)] s . We could also assume 5 0 satisfying 


f X* (E 0 r,+) H' ( Zl ) dzx 

Jr 


(n) + Sd 2 (n) + O (r 1 “), 


for certain ieK and a > 0. 

Secondly we show the solution 5 0 is in some sense unique. For this purpose, 
let us assume 5 q is another function solves CE' 0 = 0 and 


f X* (5 0 77 + ) H’ (zx) dzx = dx (rx) + S'd 2 (n) + O (rq Q ) 
Jr 


for certain S' £ R and a > 0. Then the function g := — so solves Cg = 0 and 

[ X* (gg+) H' ( Zl ) dzx = (S' - S) d 2 (rx) + O ( r~ a ) . (52) 

Jr 
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We claim that g = 0. Indeed, writing g as rj + £ (ri) H' (^i)+[?? + ^ (ri) H' (zi)] s 
f, where 

f X* ( g + f ) H' ( 2 : 1 ) dz\ = 0 

Jr 

and f = O ( r ~ a ) for some a > 0. We have 

( = Jr 'i'] + Cfdzi , 

1 (£<^ = -i^. 

Recall that in the Fermi coordinate with respect to fk,b , 

[ X* ( V + H[Cf) dz x 


A~ l d 2 I _ 1 ^ 

^ + 2 >1 dzi 2 A 2 dri 


= ~ f X* {r,+H[) 

JR 

— f X* (rj+H'i) r_1 [d ri <P*d r ri + d Zl (p*d r Zi\ dz\ 

Jr 

+ [ [-X* (rj+H[) d 2 z y + X* (ri + H[ (Su 2 - l) <p)\ dz 
Jr 


f*dz\ 


Hence f R X* ( g + H[Cf ) dzi = O (r 2 Q ) . By (l52l) , we could write £ = / 3$2 (ri)+ 
p , where /3 £ R and p (•) = O (rf “) . By the mapping property of the operator 
A'i(note that fli is in the kernel of Ad), for some cr £ (0,1), 


l/^l + 11(1 + r i) pWc 2 ’” — c 


(l +rl +a ) f g+H'Cfdz! 

Jr 


CO,- 


(53) 


On the other hand, by the equation (Cf)' L = —K 2 £, 


Cf=-K 2 £ + (Cf) 11 . 

Note that due to the orthogonality condition, {Cf)^ is small compared to f. 
Hence by the a priori estimate of the operator £, suitable weighted norm of f 
could be controlled by 


0 (|/ 3 |) + 0 11(1 + r i) pWc 2 -” • 

This together with (l53l) yields that g = 0, which implies Sq = So- Hence the 
operator £ : Si ® V —> S 2 has at most one dimensional kernel. 

Step 3. To finish the proof, it remains to show that the solutions H in step 
1 and So in step 2 indeed belong to the space S\ © D. 

Recall that the deficiency space V is spanned by 71 and 72 . Consider the 
function £ 71 , i = 1,2. We know that £72 £ S 2 . Hence by Step 1, one could find 
a solutions gi satisfy 

Cgi = C^i. 
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Note that 


Qi C% G 5i 


(54) 


for some constants q, where the function S is a non-decaying term introduced 
in (I5U1) . 

Hence by Step 2 and the asymptotic behavior of gi and 7 i, gi — Ji = dfBo,i = 
1, 2, for some constants d n . Then by (l54l) , 


CiXZ Ti — did. 0) * — 1,2. 


It follows that S — fcpiyi — fci ,272 £ Si and So — /e 2 ,i 7 i + ^ 2,272 £ Si, for some 
constants = 1, 2. The proof is completed. ■ 

Having proved the Fredholm property, we proceed to show that the operators 
involved are real analytic. 

Lemma 21 The map N is real analytic. 

Proof. This follows from the fact that A is a linear operator and the func¬ 
tion u 3 — u is a real analytic function. Note that the subtle point here is that 
the definition of N involves the diffeomorphisms 4> Sj t, Sj t. These maps are cer¬ 
tainly not real analytic with respect to the r, z variables, since there is a cutoff 
function appeared in their definition. However, these diffeomorphisms are in¬ 
deed real analytic with respect to the parameters s and t, which could be seem 
from the explicit expression m of the Fermi coordinate(Notice that / depends 
analytically on k and b ). Indeed, for a point p = (r, z), by definition, 


(p) = i)p s ,t + (1 - r])p. 


Recall that the Fermi coordinate of p S: t with respect to fk+ s ,b+t is equal to 
(ri,zi), which is the Fermi coordinate of p with respect to the curve fk+ s ,b+t- 
Hence we have the relations 



and 



The real analyticity follows from these relations. ■ 

Definition 22 A solution u is nondegenerate, if and only if the linearized op¬ 
erator C : Si © V —► S 2 is surjective. 


By the results of [2] and 12!, nondegenerate two-end solutions do exist. 
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Proposition 23 The set of solutions to © satisfying (H5l) has a structure of 
real analytic variety of formal dimension 1. Furthermore, if a solution u is 
nondegenerate, then locally around u, the solution set is a one dimensional real 
analytic manifold. 


Proof. The function u s ^t,^ is a solution of the Allen-Cahn equation, if and only 
if 


N ( s, t, ip) = 0. 


(55) 


Since N (0, 0,0) = 0, we write equation (l55l) in the form 


DN (0,0,0) (s, t,ip) + f [DN (Is, It, lip) - DN (0,0,0)] (s, t, ip) dl = 0. 

J o 

Then the result follows from the fact that N is real analytic and of Fredholm 
index 1 and the real analytic implicit function theorem(for example, see [5]). ■ 


4 Analysis of solutions on the boundary of the 
moduli space 

Two different types of two-end solutions to equation © have been constructed 
using Lyapunov-Schmidt reduction method in [2] and [12l . Let us briefly de¬ 
scribe these solutions. The first type of solutions is constructed in [2] and has 
the property that their nodal curves are close to suitable scaling of a solution to 
the Toda system. We call them Toda type solutions. The growth rate of these 
solutions is close to y/2 (but greater than y/2). The second class of solutions are 
those constructed in m- Their nodal sets are close to the catenoids, which we 
know are described by the function er = cosh (ez) , where e is a small parameter. 
We call them catenoid type solution. The growth rate of these solutions are of 
the order e -1 , hence tends to infinity as e —► 0. 

As we discussed in Section 1, we expect that the moduli space of two-end 
solutions is diffeomorphic to R. In particular, we expect that there exists a 
one-parameter family of solutions, at one end of this family(the “boundary” of 
the moduli space), the solutions should be the Toda type solutions, while on 
the other end of the moduli space, the solutions should be the catenoid type 
solutions. 

In this section, we would like to analyze the solutions near the boundary 
of the moduli space. Our purpose is to prove that if V u (recall that V u is the 
intersection point of the nodal set of u with the coordinate axes) is on the z 
axis and \V U \ is large, then the growth rate of u is close to \/2. (with additional 
efforts, one could also show that u is actually a Toda type solution, but the 
proof of Theorem Q] don’t need this fact.) We shall also show that if V u is on 
the r axis and \P U \ is large then u is a catenoid type solution. 
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4.1 Analysis of Toda type solutions 

We shall first analyze the solutions whose nodal set has two components which 
are very far away from each other. We expect that these solutions are Toda 
type. Our main result here is 

Proposition 24 Let u be two-end solution. Suppose V u is on the z axis and 
\V U \ is large. Then the growth rate of u is close to \[2. 

We use q (•) = q e (•) to denote the solution of the following Toda equation: 

c 0 q” + — q' - Cl e- 2 ^ = 0, q' (0) = 0. (56) 

r 

Observe that explicitly, 


Qe (r) 


q {er) - ^lne 


for some e > 0, where q(r) = In w ith a = 2x ^ Cl . In particular, 

q e (r) — hi r—C e tends to 0 as r tends to infinity with C e a constant depending 
on e. In the sequel, we choose e such that q (0) = \V U \. 

We shall follow similar notations as that of Section [2] For example, the 
nodal curve of u in the upper half plane will be the graph of function /. The 
Fermi coordinate with respect to the graph of / will be denoted by (n, z\). We 
also have the cutoff functions ??, rj + 7 and the solution u will be written the form 
u = u + <p, where u is an approximate solution: 


u — Ti\ -F H\ T 1, 

with the function Hi defined similarly as that of Section [2] using the cutoff 
function 77 and the heteroclinic solution H. 

The main idea of the proof is to compare / with the solution q e of 
by analyzing the equation satisfied by /. The main step will be establishing 
suitable decay estimate for the function </>, as we have already done in Section [2] 
Our starting point is the fact that if V u is on the 0 axis and \V U \ is large, then 
ll/ , |li“(( 0 ,+oo)) small. This follows from an application of the balancing 

formula. 

We shall get a L°° estimate for </>. In the following, H-H^ stands for the 
L°° ((0 + 00 )) norm. To simplify the notations, we only consider the case that 
the radius of the Fermi coordinate is large enough. 


Lemma 25 Suppose u satisfies the assumption of Proposition \2f\ Then 


Moo<c 


11 

r 2 


+ Ce -W2/(°). 
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Proof. We only sketch the proof, since many computations will be similar to 
that of Section [2] 

Recall that 0 satisfies the equation 

L0 = [E(u)} 11 + [E^)] 1 + P(0). (57) 


P (c t >) is a higher order term of 0, and [E (u)]^ = E (zt)f + E (zi)| is also small 
compared to 0. As in Section [2l 


[E (u)] 1 - = [E {'HI)] 1 ' + O (.f " 2 ) + O {r~ 2 f 2 ) + O (ft" f") + O (ft'/-" 1 ) (58) 

+ O { ft " 2 ) + O {ti 2 ) + O (tif") + O ( tif") + O (e'^ 215 ) . 


Project equation (l57l) onto rfH [, we could show 


c 0 f 


+ 


C oh" 

f/2 


C 0 /' 


- (1 + o(l))cie 


-V2D 


(1 + f 2 ) 2 1 + f 2 riy /1 + / ,: 

= O { ft "/") + O {h' 2 ) + O {tif") + O (. tir _1 ) + O {tif") 


(59) 


+ o (|| 0 * (n, OIlL) + O (110* (n, •)IL + o ( 110 * (n, OIL 11^0* (n, OIL) 

+ o(||0 ri 0* (ri ,.)|| oo e-^))+o(||0 I . 1 0* (n.OIlL) +0(110* in, OIL IK0* (n.OlIJ 
+ O if" 11^0* (n, OIL) + ° (/" 11^0* (n, OIL) + 0 (/'" 11^0* (n, OIL) 

+ O (r- 1 ||0 ri 0* (n, OIL) + ° ( r_1 /' 113* 0* (ri, OIL) ■ 


Insert this into (1551) , we obtain 


(^(S)) X 


= [P {HI)] 1 + O (r" 2 /' 2 ) + O {h'r- 1 ) 

+ O {ti 2 ) + O {tif") + O (e~' /2D ) 

+ O (ft'ft") + O {hf") + O (ft" 2 ) 

+ o(||0* (n,OllL)+0(11^0* (n.OIlL) 
+ 0(110*0* (n.OIlL) ■ 


Then the a priori estimate of the operator L tells us that 


II0IL < o 


11 

r 2 


+ Ce - 2 V2f(0)_ 


Our next aim is to estimate the L°° norm of 


f(r 


Lemma 26 Let e be introduced above, then 


f 


+ ||/"|loo < Oe 2 . 
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Proof. The starting point of the proof is still equation (15U1) . Applying Lemma 
l25l we infer from (1551) that 


C 0/" , C 0f' _ c V2D 

(1 + /' 2 ) 1 n(l + /' 2 )^ 

= o( £ 2 ). 


(60) 


Let to be a point where 


to 



f' 

This point exists because ^ —> 0, as r —> +oo. 
At the point to, by (l6()l) , / will satisfy 


co/"(to) + co^-c ie -^°) 
to 

= O (e 2 ) • 


(61) 


We claim 

f" (to) > 0 . (62) 

Indeed, if this is not true, then due to the fact that /" (0) > 0, there will be 
another point t\ < to such that /" (ti) = 0, and /" ( t ) < 0, t £ (ti,to) • Since 
/' > 0 , we find 


/'MV 


< o, r e (ti,t 0 ). 


This contradicts with the fact that t-&l 

to 

be concluded that 


L 


. From (RJT1) and (KHZll . it may 

OO 


fM = 2l e -V2D(t 0 ) _ f n (t()) + Q / 2) 
to Co 

< Ce- 2v/2/(0) < Ce 2 . (63) 


Here we have used the fact that D (to) > 2/ (0). This proves estimate for the 
L°° norm of t -. The estimate for /" is a direct consequence of (1M1) and (1CT1) . ■ 
Let p = f + \J 1 + f r2 h. Next we show that p is indeed close to the solution 
q of the Toda equation in a large interval. Set b = b e = | ^ | . 

Lemma 27 There exists a > 0, such that 


I P (r) -q(r) | < Ce a , r <E (0,6). 

Proof. Denote by e the function p (e _ 1 r) + lne — q (r). Then e (0) = O (e 2 ) 
and e' (0) = 0. By the previous lemmas, < e 2 . Using this fact, we see that 

the function p = y/l + f' 2 h, + f satisfies the equation 

CoP " 3 +-_ Cl e~^ D = O (e 2+ “) . 

(1+p' 2 ) 2 ri(l+p' 2 ) 2 
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In particular, this combined with the fact that D — 2/ = O (s a ) yields 

p" + — - e~ 2V2p = O (e 2+a ). 
r i v ' 

From this equation, we infer that in the region where e is o (1), e satisfies 


e" + e~ 2V2<l e = O (e 2 ) + O (e a ) . 

The conclusion of the lemma then follows from the variation of parameters 
formula. ■ 

By definition q (r) = q (er) — lner, hence rq 1 ( r ) = erq' (er), which implies 
that 

bq ( b) = V2 + o ( 1 ) . 

Notice that in Lemma [271 we actually could also estimate |e'| < Ce a . From this, 
we infer that 

bf(b) = V2 + o{ 1 ). (64) 

Now we are in a position to prove the main result of this section. 

Proof of Proposition 1241 Similar arguments as before yields 

\4>{r,z)\<C-±-s+Ce- D W. 

Equation (15U1) then becomes 

-+---- = (1 + o ( 1 )) e ~^ D ^ + O (rq 3 

( 1 +p ' 2 ) 2 ri(l+p ' 2 ) 2 


Integrating from t 0 to t\ leads to 

tlP ' {tl) , - top,{to) t = (! + o (1)) T e-^sds + O (t^) . 

(l+p'Ctr ) 2 ) 2 (l+p'(t 0 ) 2 ) 2 Jt ° 

This together with (1SH) tells us that 

p' ( r ) r > \/2 + o (1) for r > b. 

Let 6 > 0 be a fixed small constant. We claim that when e is small, 

p' (r) r < a/2 + 6 for r > b. (65) 

Indeed, suppose ( 6 , b*) is the maximal interval where p' (r) r < \[2 + S. Then in 
this interval, elementary geometrical facts implies that 


D(r) > v / 2|lne| + 2q(6) + 2 
= V2 |In + 2 q ( 6 ) + 2 




s 


(In r — In b) — C. 
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Therefore we have 


, _ffl-r = (1+0(1)) r + 0 ((,-') 

(l + p'(&*) 2 ) 2 (l I-i V) ; 

= 0(b-‘). 

This implies that when e is small, b* = +oo. Applying (1551) , we finally get 
lim r _s.+oo p' (r) r = y/2 + o (1). The proof is thus completed. ■ 

4.2 Uniqueness of catenoid type solutions 

Beside the planes, catenoid is the first example of embedded minimal surfaces 
with finite total curvature; it is rotationally symmetric with respect to its axis 
and the only minimal surface of revolution(up to a homothety). In the (r,z) 
coordinate, the one-parameter family of catenoids C E could be represented by 
the function 

er = cosh (sz ), 

with e > 0 being the parameter. As we mentioned before, in |12j . for each e 
sufficiently small, a solution u e of the Allen-Cahn equation is constructed. The 
nodal set of u e is close to the catenoid C e . In particular, V Ue is on the r axis, 
\V Ue \ is of the order O (e _1 ) , and the growth rate of u e is also of the order 
0(^). 

In this section, we wish to prove that this one-parameter family of solutions 
u e is unique. This is the content of the following 

Proposition 28 Let u be a two-end solution of the Allen-Cahn equation. Sup¬ 
pose V u is on the r axis and \V U \ is large. Then there exists a small e' > 0 such 
that u = u E >. 

The nodal curve Af u in E + could be written as 

{{r,z) :z=f(r)}, 


where / is a function 

/ : [I’Pul, +oo) —> R. 

Note that by the results in the previous section, / is asymptotic to ci In r + C 2 
as r —> Too for some constants ci and c 2 . One could also write Af u D E + = 
{(r,z) : r = g (z)} for some function g : [ 0 , +00) — > R. 

For convenience we introduce the parameter e = \V U \ 1 • Observe that by 
the assumption of Proposition [28l e is small, thus by the validity of De Giorgi 
conjecture in R , locally around the nodal curve, u looks like the heteroclinic 
solution. 

Let l be a large but fixed constant. As a preliminary step, we would like to 
get some rough information about the slope of the function /. 


44 




Lemma 29 We have 


f (r) > j +o(l), r G [e 1 ,le x ] . 


Proof. Let X = (0, 0,1) be a constant vector field. We have 

rle~ 

Jo 

On the other hand, 


\u 2 r + F ( u ) 


rdr > Cle . 


/ { 

]^u 2 r +F{u) 

/ann{z>o} l 



X ■ v - (Vu • X) {Xu -v)\dS 
f (is -1 ) Is -1 + o ( 1 ) le -1 . 


Combine these two estimates and use the balancing formula, we get the desired 
result. ■ 

To get more precise information, we again need to work in the Fermi coor¬ 
dinate (s,t) around the nodal curve A f u , where s is the signed distance to A f u 
and t is a parametrization of A f u . 

We slightly abuse the notation and still use B u to denote the maximal domain 
where the Fermi coordinate is well-defined. Similarly as before, let g be a 
cutoff function supported in B u , with Xg supported near the boundary of B u . 
Introduce the approximate solution 


u (r, z) 


rjH + (1 — 77) 


H 


where Ti is defined through X*'H{s,t) = H{s — h{t)). Here X is the map 
(s, t) —> (r, z). Write u = u+(f>. The small function h is chosen such that </> 
satisfies the orthogonal condition 


X* (rf<j>W) ds = 0, 


where X*^' (s, t) = H' (s — h (t)). 

To analyze the solution u, we firstly study A f u in the region where r G 
(e _ 1 ,l£ _1 ) . In this region, many calculations will be explicit to in the Fermi 
coordinate (x,y) with respect to the graph of the function g, which is 

/ ' = + 

l Z = V ~7^- 

Hence x is the signed distance and y is a parametrization of the curve. 

We shall estimate the L°° norm of the perturbation term (j). 
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Lemma 30 Suppose u satisfies the assumption of Proposition I SSI and </> is de¬ 
fined above. Then 


Moo < Ce 2 . 


Proof. We only consider the case that the size of the Fermi coordinate is large 
enough. The general case could be handled using the arguments of Section [2] 
Again many computations here are similar as before. 

We need to analyze E (u). First of all, consider the case r £ (e^ 1 , /e _1 ) . In 
the Fermi coordinate, the error of the approximate solution has the form 


E(u) 


H"h' 2 - H'h" 

A 


d x A 
2 A 


+ 


dyA 
2 A 2 



H'h'. 



H' 


Keep in mind that H' is evaluated at x — h {y ), not x. By Lemma [?H \g'\ < C, 
hence 


d r x 

r 


l 

y/l+g ' 2 

9 + 7 ^ 


1 

gV 1 + g ' 2 


X 

g 2 (1 + g' 2 ) 


+ o Or 3 ). 


From | g"\ = o(l), we infer that 


d x A 
2 A 


(i + g ,2 ) § 


W'fx 

(i + g' 2 f 


+ 0 ( g " 3 ). 


It follows from these expansions that the projection of E (it) onto TL' has the 
form 


[ X* (r)TL'E (■u )) dx = — 

Jr 


cq g 


c 0 


-h" 


(i + g ,2 Y 
H’ 2 
~A + 


v/l + 


g 

O ( h ,2 ) 


5 

0(h'g") 


+ O ( h'g'") + O ( g " 3 ) + O (g~ 3 ) 
+ O (. hg ,/2 ) + O (hg 2 ) . 
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Set h {y ) = y/l + g' 2 h (y) and p ± (y) =g(y) + h (y ). Then 
g" 1 h" 

1 + g ,2 + g y/i + g ' 2 

g" ^ 1 h" 2 ti ( 1 Y 

i + g ,2 + g i + s ' 2 ^i + 5 ' 2 Vv/i + ff'V 

* ( * y 

= — y 9 H-hO (g"h) + O (. g"h') + O {g"'h) + O (, g"'h') 

+ O (. g"h ) + O ( 5 "h') + O (h 2 ) + O (/i/i') + O (g^h). 

It follows that 

\/l + 3 /2 f X * {rfH'E (S)) dx = --^4 + - 
Jr 1 + Pi Pi 

+ O (<A) + O (. g"h ') + O {g"’h) + O {g'"ti) 

+ O {g"h) + O ( g"h ') + O (h 2 ) + O (hh') + O {g~ l ti) 
+ O {ti 2 ) + O ( g " 3 ) + O ( g ~ 3 ) + O {h'h"). 

On the other hand, we could also show that 

[ X* (r(H'E{u))dx = o(||0* (y, -)\\ C i)+o (\\d y <j>* (y, •)|| c . 0 )+o (||S 2 0* (y,-)|| c0 ) • 

(66) 

Here the norm is taken as a function of x variable. As a consequence of the 
above two equations, we have 

-TT~n +~ = 0 {ti 2 ) + O {tig") + O {tig"') + O {g " 3 ) 

+ O {g - 3 ) + O ( hg" 2 ) + O ( hg ~ 2 ) + O {h'h") 

+ 0 ( 11 ^* ( 2 /,-)llcO + °{\\ d v ( t ) * (y.Ollco) + o(|j<9 2 0* (Pr)|| C o) • 

(67) 

Next we consider the case r > e _1 L In this case, one could use the Fermi 
coordinate (ri,Zi) with respect the curve z = f (r). The analysis of E {u) in 
this case is almost same as that of the previous sections and we omit the details. 
Now we write the equation satisfied by 0 into the form 

LJ>=[E{u)] ± + [E(u)}"+P (<(>), 

where 

[E («)]" = X j ^ ^^ dX gW and [E (u)]- 1 - = E (u) - [E (u )] 11 . 
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In terms of Fermi coordinate, in the region where \g'\ < C, 

[E(u)] ± = o(\p f ;\ 2 )+o(pi 2 ). 


In the region where \ f'\ < C, [E (m)] ± = O (rq 2 ) . By the a priori estimate of 
L , we could obtain 


Moc< [E(u)] 


< Ce 2 . 


We emphasize here that 


[£(«)] X 


should be estimated in the whole plane. It 


is also worth mentioning that the term O (\p '{| 2 ^ in [E (it)] 1 " should be handled 


using equation (El). ■ 

Our next aim is to show that in the interval (e _ 1 ,Zgr _1 ) , the function r = 
Pi(z) is close to the function r = e _1 cosh (£ 2 ). This is the content of the 
following 


Lemma 31 For z £ (0, / (7) £ x ) , 


\p\{z)—E 1 cosh (ez )I < Ce. 


Proof. From equation (l67l) and the estimate of <f>, we deduce that the function 
Pi satisfies 


p" 1 
1 + pi 2 pi 


O (£ 3 ) . 


At this stage, we introduce the scaled function pi (z) = epi (£ 1 z) . Then 


p’ 1 (z)=p’ 1 (e^z), 

p'/(z)=£-V/(£-^). 


It follows that 


Pi - = o (£ 2 ) , Z€(0,l). 

Pi 


Observe that the function cosh z satisfies the equation 


(cosh^) ,/ — 


1 + (cosh z) 
cosh z 


/2 


= 0 . 


Let u) (z) = pi (z) — coshz. Then 


( 68 ) 


(69) 


w ( 0 ) = pi ( 0 ) - 1 = epi ( 0 ) - 1 

= £ (e^ 1 + h (0)) —1 = 0 (e 2 ) . 

We claim that |w( 2 )| < Ce 2 for 2 £ (0,1). Indeed, subtracting equation (1551) 
with (l69l) , we get 

w" — 2 tanh zcu' + uj = O (e 2 ) + O (w 2 ) + O (u/ 2 ) := ifi ( z ). 
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Let £1 and £2 be two linearly independent solutions of the homogeneous equa¬ 
tion: 

£■' - 2 tanh + & = 0, i = 1,2. 

Explicitly, we can choose 


£1 ( 2 ) = cosh' 2 : = sinh 2 , 

£2 (~) = d e (e” 1 cosher) | e= i = — cosh 2 + 2 sinli 2 . 
The Wronsky of these two solutions are 


W(z) := 


sinh 2 cosh 2 

cosh 2 + 2 sinh 2 2 cosh 2 


= cosh 2 2 . 


By the variation of parameters formula, 


“ (z)= 6 w L w{.) ds ~ (iiz) l — m~ +0(ey 


The desired estimate comes from this formula. ■ 

Proof of Proposition l28l Let P 2 = f + \/l + f r2 h. For r > ie _1 , projecting 
E (u) on rfH' and perform similar calculation as in Section [21 we could estimate 
the perturbation 1 in algebraically weighted norm (Remember that h could be 
controlled by 0). This leads to the equation: 


( r lP2 \ 

[V^TW) 


= o 


1 + e 2 r 2 


Introduce the scaling of p 2 : 


P2 (o) = eP2 (e V) ■ 


We find that in (/, + 00 ), p 2 satisfies has the equation 


( r ip' 2 \ 

V \A + P2 j 


= 0 



(70) 


From this equation, we deduce 

lim np' 2 (r 1 )-lp' 2 (l) = 0(e 2 ) . 

7*1—>-+00 

On the other hand, by Lemma [dll 

lp' 2 ( l)-l = 0 (e 2 ) . 


Hence 

lim rip' 2 (ri) = 1 + O (e 2 ) . 

7*1—>-+00 V 7 
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Consequently, the growth rate of u is equal to 


lim np' 2 (ri) = - + O (e) := 

Obviously, e' = e + O (e 3 ) . 

Now let r = e ,_1 cosh (e'z) be the catenoid with the same growth rate as 
u. Then there is a solution u e ' with the slope e ,_1 whose nodal line is close 
to this catenoid(the nodal line in E + decaying algebraically to a vertical small 
translation of the catenoid. Actually, the error could be estimated by some 
positive power of e. Our aim is to show that u = u E i. To see this, we need to 
estimate the function pi (z) — £ ,_1 cosh (e'z) and P 2 (ri) — e' -1 cosh -1 (e'n). 
Note that 

ee' -1 cosh (e'e~ 1 z) — coshz = (l + O (e 2 )) cosh (l + O (e 2 ) z) — coshz 

= 0 (£ 2 ). 

Hence using Lemma m we obtain p\ ( z ) — ee' 1 cosh (e'e 1 z) = O (e 2 ) . This 
implies that 

Pi (z) — £ /—1 cosh (e'z) = O (e) . 

Using COD , we then find that P 2 (?r) — e' 1 cosh 1 (eVi) = O (e). 

Therefore the nodal line of u is close to the nodal line of u e > at the order O (e) . 
Then, applying the mapping property of the Jacobi operator of the catenoid, a 
contraction mapping argument shows that u = u e '. We refer to m for similar 
arguments in the case of four-end solutions in R 2 . ■ 


5 Concluding the Proof of Theorem |1] 

In this section, combining the results of the previous sections, we would like to 
finish the proof of Theorem [I] 

Let M be the set of all two-end solutions to the Allen-Cahn equation with 
growth rate larger than i/2. Consider a catenoid type solution u e arising from a 
largely dilated catenoid. As we mentioned in Section [3l u e is nondegenerate.By 
Proposition [231 locally around u e , the set of two-end solutions is a one dimen¬ 
sional real analytic manifold, which we denote as the image of a map g 

6 '■ ( — £> $) M, 

for some small S > 0. Using the compactness result in Section [2] and by the 
structure theorem (1191) . g has a global continuation: 

g : (-J, + 00 ) —> M. 

We claim that the growth rate of the solution g(t) tends to \/2 as t —> + 00 . 
Indeed, as t — > + 00 , V e ( t ) could not remain bounded. (Recall that V g ( t ) is 
the intersection of the nodal set of g(t) will the r or z axis.) Otherwise by 


50 


the compactness, the image of g will be a closed loop, which could not be 
true since the family of catenoid type solutions u e are not compact. Also, as 
t —> +oo, V e (t) could not be on the r axis, this follows from the uniqueness 
of catenoid type solutions, Proposition [5S] Therefore, P e (t) will be on the z 
axis and |’P e (t)| —> +oo, as t —> +oo. By the analysis of Toda type solutions, 
Proposition [2H the growth rate of the solution g(t) will go to \/2 as t —> +oo. 
This finishes the proof. 


6 Appendix 

6.1 Monotonicity of two-end solutions 

In this appendix, we sketch the proof of monotonicity of two-end solutions(Proposition 
[3]in Section[2]) using the moving plane argument. Essentially, we follow the proof 
of monotonicity for four-end solutions of 2D Allen-Cahn equation in [13] , But 
their is a slight difference here. Namely for the two-end solutions in dimen¬ 
sion three, to start the moving procedure at infinity, one need to have suitable 
control of the asymptotic behavior of the solution(estimate (1731) below), while 
this is not needed in dimension two case. The asymptotic expansion we need is 
provided by the results of Section [2] 

Proof of Proposition [3} Let it be a two-end solution. We first prove its mono¬ 
tonicity in the r direction. To use the moving plane machinery, we will work in 
the usual Euclidean coordinate (x. y , z) . Set U (x. y , z) := u (^/x 2 + y 2 , z'j = 
u (r, z ). 

Suppose the growth rate of u is equal to k > \[2. Hence the asymptotic 
curve of its nodal line in the first quadrant is 

z = k In r + b 


for some b € R. 

Let £o > 0 be a parameter and define 

U {x, y,z) = U (x, y, z\ x 0 ) ■= U {2x 0 - x, y, z). 

Certainly U (xo ,y,z) = U (xo,y,z ). Note that U actually depends on the pa¬ 
rameter Xq. 

The first step is to show that the moving plane procedure could be started 
at + 00 . We claim that for xo large enough, say Xo > a o, 

U (x, y, z ; x 0 ) < U (x, y, z ), for x < x 0 . 

First of all, we consider the region where \U\ is not close to 1. Recall that the 
nodal set of U is close to 


z = k In (2xo — x) 2 + y 2 + b. 
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We have 


In \J ( 2 xq — x ) 2 + y 2 — In \Jx 2 + y 2 
4®q — 4xqX \ 


= iln|l + 




(71) 


Note that for x < Xq, 4j ° j x ° x > o. Fix a small positive constant e. 

If > e, then by O , 

In \J^2xo^x)^~-\-y 2 — In \J x 2 + y 2 > ^ In (1 + e). 

By the results of Section [2] for r large, we have 

u (r, z)-H (zi) = O (r -2 ) . 

Here z\ is the signed distance of (r, z) to the nodal line. As a consequence 
U (x,y,z;x o) - U ( x,y,z ) 


(72) 

(73) 


= H 
- H 


( 2 xo - xf + y 2 , z'j -u (r, z) 

^ z — k In \J (2®o — x ) 2 +y 2 — b cosdi^j 
^ z — k In \Jx 1 + y 2 — b cos 62 ) + O ( r~ 2 ) , 


where = O (r 1 ) . It follows from (l72l) that 


U (x, y, z; xq) - U (x, y, z) < H - k In \J (2x 0 — x) +y 2 -b 

— H (z — k In \Jx 2 + y 2 — bj + O (r~ 2 ) 
< 0 , 


for r large enough. 

If 4^-4xox e(0j£))then 

In \j (2xo — x) 2 + y 2 — In ^/^Ty 2 = ^ ~^ X + O ^0 fro- »o) 2 ^ _ (74) 

There are two possible cases. Case 1: x £ (—oo,cco — 1) • In this case, by (17H) , 
In \J (2x 0 - x) 2 + y 2 - In \Jx 2 +y 2 > 

This together with ( 1731 ) implies that for Xo large, U < U. Case 2: x £ (xq — 1 , xo) . 
In this case, by the estimate 

d r (u (r, z)-H ( 21 )) = O ( r~ 2 ) , 
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we find that (here one should also use some estimates of the Fermi coordinate) 
U — U < -C x2 ° ~ 2 XoX + O (r- 2 (xo - X)) 

for certain positive constant C. Hence if we choose Xo large, U < U. Therefore, 
to prove the claim, it remains to consider the region where \U\ ~ 1. Observe 
that for Xq large enough, in the region where U ~ 1, U is also close to 1. Now 
let <f := U — U. Then in the region x < Xq, f satisfies 

~^(x,y,z) i p + (u 2 + UU + U 2 — l) f = 0. 

Since limsupuj.^ ^.)i^. +00 f < 0, by the maximum principle, tp(x,y,z) < 0, for 
x < Xq. This proves the claim. 

In the second step, we define 

x* = inf {x : U (x, y, z ; xq) < U (x, y, z ) for x < Xq and Xo £ (x, ao)} . 

We show that x* =0. To see this, we first prove that r* < ao. Indeed, by the 
first step, 

U (x, y, z; a 0 ) < U {x, y , z) for x < a 0 , 

U ( a 0 ,y,z;a 0 ) = U ( a 0 ,y,z ). 

Hence by the Hopf Lemma, d x (U (•; ao) ~ U) > 0 for x = ao. Then standard ar¬ 
guments together with the asymptotic behavior of U implies that U (x, y, z; xq) < 
U, for x < Xo and xq sufficiently close to ao- Then one could use this type of 
arguments to show that x* = 0. (The plane could be moved to the left until the 
inequality U < U is violated at infinity.) 

The monotonicity in z direction could be proved similarly using moving 
plane. This completes the proof. ■ 
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